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ABSTRACT 



Aims. We present the results from high-resolution spectroscopic measurements to determine metallicities and activities of bright stars 
in the southern hemisphere. 

Methods. We have measured the iron abundances ([Fe/H]) and chromospheric emission indices (logi?[jj^) of 353 solar-type stars with 
V=l. 5-9.5. [Fe/H] abundances are determined using a custom ;^-^ fitting procedure within a large grid of Kurucz model atmospheres. 
The chromospheric activities were determined by measuring the amount of emission in the cores of the strong Can HK lines. 
Results. The sample of metallicities has been compared to other [Fe/H] determinations and was found to agree with these at the 
±0.05 dex level for spectroscopic values and at the ±0. 1 dex level for photometric values. The distribution of chromospheric activities 
is found to be described by a bimodal distribution, agreeing well with the conclusions from other works. Also an analysis of Maunder 
Minimum status was attempted and it was found that 6+4 stars in the sample could be in a Maunder minimum phase of their evolution 
and hence the Sun should only spend a few per cent of its main sequence lifetime in Maunder Minimum. 

Key words. Stars: abundances - Stars: activity - Stars: atmospheres - Stars: planetary systems 



1. Introduction 

High resolution spectroscopy is an important tool used in the 
understanding of stars in the local neighbourhood and can 
give an insight into various fundamental properties of the 
stars themselves (e.g. lEdvardssonet aPllQQSt IWrightl l2004aL 
IValenti & Fischeij l2005l) . Indeed, knowledge of stellar funda- 
mental properties can also have a direct bearing on other facets 
of stellar astrophysics such as the potential for discovering 
planetar y systems. After only three exoplanets had been dis- 
covered JGonzaiea (Il997h noticed all were orbiting stars with 
a higher metal content than the Sun. This trend has contin- 
ued and been confirmed by various authors, most recently by 
iFischer & Valentil (l2005h who have shown that the probability of 
discovering an exoplanet around a solar-type star is proportional 
to the stellar metal abundance following a power law described 
by 0.03x1 O^"'^*^'^', consistent with the predictions made by the 
core accretion models of planet formation. Along with a star's 
metallicity another important feature is the star's chromospheric 
activity. 

Stellar chromospheric activity is intimately r elated to the 
stellar dynamo, magnetism and the stellar rotation (Middelkoon 
ll982al: lMiddelkoop 1982b; Rutten 1984) . Since the measure"- 
ment of activity gives us a feel for the motions on both large 
and small scales in a stars atmosphere, it can be used to se- 
lect the most inactive and hence stable stars to scrutinize for 
exoplanets. It has been found that chromospheric activity can 
mask or i n some cases exhibit planetary radial-v elocity signa- 
tures (e.g. lOueloz et alj|200ll: [Henry et al.ll2002l) . Therefore it 
seems prudent for any radial- velocity planet search program that 
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is aiming to reach precisions of a few ms"' to select the most 
inactive stars. We have begun a radial-velocity project to de- 
tect both hot Jupiter-type exoplanets and lower-mass hot Saturns 
and Neptunes that have the potential to transit their host star. In 
Section 2 we discuss all observations and reduction methods. In 
Section 3 the methodology for extracting the spectroscopic in- 
formation is explained: both the activity analysis and the stellar 
photospheric properties such as metallicity abundance. Section 4 
compares the values determined here with those of other works 
in the literature. In Section 5 the results are discussed in terms of 
both the current status of exoplanets and stellar evolution. Lastly, 
in Section 6 we discuss the conclusions of the project and plans 
for future work using this dataset. 

2. Observations & Reduction 

All target stars and calibration data were observed using 
the Fibre-fed Extended Range Optical Spectrograph (FEROS) 
mounted on the MPG/ESO - 2.2m telescope on the La Silla 
site in Chile. To gain access to the entire southern sky the 
run scheduling was split across two independent observing runs 
seven months apart. The runs were three nights in length, from 
the 2"'' to the 5* of February and September 2006. This cul- 
minated in over 350 stellar observations with exposure times in 
the range 120-480 seconds. These exposure times gave rise to 
S/N ratios of 100-200 in the continuum around the iron line at 
7500A and a median of -60 at the CaHK fines (3955A), witfi 
only a small tail of low S/N stars reaching down to values of 30. 
All calibration files needed for the reduction of the stellar spectra 
(flat-fields, bias and arc frames) were obtained at the beginning 
and end of each nights observing, following the standard ESO 
calibration plan. 

The reduction of all the spectra followed the standard reduc- 
tion techniques. Firstly the bias signal was removed from each 
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Fig.l. The extracted FEROS spectra for the star HD59100 
around the Can K and H hnes. The star has a spectral type of 
G2V and is chromospherically quiet (log/?[|j^=-4.92). All the 
four passband regions (V,K,H&R) have been highlighted on the 
plot. The spectrum has been normalised to the continuum region 
3991-401 lA which represents the R passband. All other features 
are relatively much weaker absorption lines, not noise. 
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Fig. 2. This plot shows the linear fit applied to the FEROS spec- 
tra to convert the dataset to the MW system of measurements. 
All stars used as calibrators were taken from the MW project. 
The fit shows an increasing linear trend with increasing S'-index. 
The overall scatter of the points is highlighted on the plot, and 
with a value of 0.0 6, this is in agre ement with the findings of 
iTinnev et all (|2002|) . IWrighj (l2004ah and J06 who all calibrated 
using high resolution data and setups similar to the MW project. 



individual frame. The overscan region was then subtracted from 
all frames individually and then trimmed off the image so as 
not to confuse any of the extraction algorithms. FEROS echelle 
data shows strong curvature in the 2D echelle image and, due 
to the use of an image slicer, has a complex slit cross-profile. 
The curved echelle profile was straightened by fitting polyno- 
mials along the orders. These were then precisely traced using 
a well exposed Halogen flatfield image to follow the lamp light 
path. Due to the half moon-like pattern of the sliced beams the 
orders were smoothed in the spatial direction to create a sin- 
gle smoothed order. This allowed the tracing algorithm to prop- 
erly follow the center of the order and not vary between the 
two sliced peaks. The traces were then clipped to remove any 
stray pixel counts and to help tighten up the trace and then the 
object and background apertures were selected. To correct for 
the CCD pixel-to-pixel response we obtained lO-i- flatfields with 
S/N ratios close to 1000 each morning and evening. Each pixel 
was median combined together to create a master flatfield frame. 
The master flats were then used to calculate the balance fac- 
tors by smoothing in the dispersion direction and dividing the 
smoothed flat by the master flat. These balance factors were 
then multiplied into each science frame to flat field each im- 
age. ThAr+Ne and ThArNe lamps were used to wavelength cal- 
ibrate all the data. The sky background was negligible; however 
the scattered light was still removed. This was done by setting 
the dekker limits large enough to include wide inter-order spac- 
ing regions and then a low-order polynomial was used to sam- 
ple any gradients along the orders. A second low-order polyno- 
mial was then used to fit to the total intensity over the whole 
order to better model the scattered light. The profile of each or- 
der was measured by sub-sampling each order profile individ- 
ually using Gaussians. This model was used to extract the ob- 
ject with an optimal extraction algorithm. The data were then 
binned to a linear wavelength step of 0.03A/pixel, or a mean 
velocity resolution of 1.42kms" '/pixel, and shifted into the rest 
frame by cross-correlating with an observation of HD 102 11 7 
(G6V), which has been shown to exhibit a radial-velocity varia- 
tion <20ms ' jJones et al.ll2002l) . 



3. Methodology 

3.1. Activities 

The methodology used to extract these high resolution ac- 
tivity indices is similar to that employed in pr evious works 
by the Anglo-Australian Planet Search e.g. iTinnev et al.l 
(2002); Jenkins et al. (2006) (hereafter J06) and Mount Wilson 
(.Duncan et al...l991.) projects. Four passbands situated in a well 
concentrated region of the UV were used for flux measurements. 
Fig. [U shows the part of the spectrum used for all stars. This 
particular spectrum is of HD59100, which is a similar star to 
the Sun with a spectral type of G2V and has a final activity of 
log/?Jjj, = -4.92. All the four passbands used to make the activ- 
ity measurements have been shown for reference. Each of the 
individual lines apparent in Fig.[T]are relatively weaker spectral 
features, not noise. 

3.1.1. Smw Calibration 



ITinnev et all (|2002|) . J06 and IWrightl ( l2004al) have shown that 
linear calibrations are needed to calibrate high re solution activ- 
ity indices onto the MW system of measurements dDuncan et al.l 
1991). The setup employs four separate passbands centered on 
and around the Can HK lines. The passbands used to estimate 
the mean flux in the continuum (V & R) are both square and cen- 
tered at wavelengths of 3891 A and 4001 A respectively. Whereas 
the passbands used to obtain the mean flux in the K and H line 
cores have triangular profiles with full-width half maximum's 
(FMHM) of 1.09A and are centered on the lines themselves at 
wavelengths of 3933. 664A and 3968.470A respectively. 

The following relation (Eq. [TJ, is used to determine the 
FEROS activity indicator S'feros: 



Nh+N^ 



K 



I FEROS 



A^R -hA^v 



(1) 



Here the A^, is the number of counts in each bandpass (where 
/=H, K, V and R) and hence this provides a measure of the ratio 
of the flux in the Ca ii line cores to that of the continuum. When 
this index was created for use at Mt. Wilson it was tailored to- 
wards use with the photometer and passband wheel there and 
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also employed a normalising factor, however by use of a similar, 
yet artificial setup, this procedure works very well for extracting 
activity indices of bright solar- type stars. 

A number of stars from iDuncan et alj ( 1199 lb were selected 
for use as calibrators onto the MW system, these are shown in 
Table. [3] For any activity study, the selected comparison stars 
should not introduce any offset into the final 5-indices, how- 
ever this remains difficult due to the lack of long-term moni- 
toring data in the literature for stars in the southern hemisphere. 
Two of our calibration stars, HD1835 and HD10700, have re- 
cently had their long-term va riability and activity trends pub- 
fished. JLockwooditlD dZQOTl) have shown that HD1835 did ex- 
hibit variation at the level of ~+0.05 between 1985 and 2002. 
However, the star only spent a limited time at these extremes 
and was most ly found to be only ~+0.02-0.03 from the mean 
activity value. iHall et all (12007) found a mean S'-index for this 
star of 0.362 using 82 observations over a period of ~11 years, 
and this mean is 0.012 larger than the final computed value here. 
HD 10700 was also measured by Hall et al. and they found the 
mean 5-index to be 0.175 over 10 years and 67 observations, 
which is close to the value of 0.179 found here (i.e. two mea- 
surements of 0.178 and 0.180) within the uncertainties. Fig. |2] 
shows the linear fit between the FEROS and the published MW 
values. The fit has a gradient of 1 . 144+0.001 and an intercept of 
0.018+0.006. The overall scatter of the points around the linear 
trend is 0.006, which is similar to the scatter found by J06. This 
low level of scatter gives confidence in the robust nature of the 
reduction and analysis procedure. The calibration error is found 
to be consistent with the removal of any of the data points. 

IWrightl(l2004al) and J06 used measurements of the stable star 
T Ceti (HD 10700) as a proxy for the random uncertainty intro- 
duced in the reductio n procedure. As t Ce ti has been shown to 
be extremely stable (iBaliunas et al.|[l995l) this is a useful tool 
for estimating the uncertainty introduced in the reduction proce- 
dure. Only two measurements of this star were possible due to 
the star's position and bad weather, and hence it could only be 
used as a weak uncertainty estimator The measurements were 
obtained over two separate nights and have shown the stability of 
the reduction. The computed S'-index for both measurements are 
0.178 and 0.180, indicating that the procedure is robust. As the 
reduction uncertainty is weakly estimated here, the RMS scatter 
can be used as a first approximation of the t otal uncertainty , and 
also by comparison with the findings from lWrightl ( l2004al) and 
J06 where both find uncertainties of ~3-5%. We estimate that 
the total error included in each individual S activity indicator is 
around 3-5%. 



3.1.2. logR^^ Fit 

The measured S'-index contains photospheric information and 
for activity analysis one wants exclusi vely chromospheric infor- 
mation, therefore the methodology of Noyes e t al.l (1 1984 a) has 
been used to remove the photospheric component of the activity 
index. This method normalizes the flux to the bolometric lumi- 
nosity of the star by use of empirical relations (see Noyes et al. 
for equations). Fig. [3]shows the final logR'^^ feros values for all 
stars used to calibrate onto the MW system of measurements. 
The plot highlig hts the close re l ations hip between the published 
MW values from lDuncan et al](ll99lb and these values. A linear 
best-fit trend is applied to the data and this is shown by the solid 
line in the plot. The gradient of this trend is 0.981+0.003 and 
has a negative intercept of -0.092+0.250. The slope of this fit is 
in agreement with that found by J06. No other calibration is nec- 
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Fig. 3. The final Xq^R'^t^ values for th e FEROS cal i bration 
dataset against the values published in i Duncan et al.l (Il99ll) 
(MW). There appears a tight fit around the 1:1 relation, which 
is shown by the dotted line. The dashed line represents the best- 
fit linear trend to this data and the tight relationship means no 
further calibration is required. The RMS scatter around the fit is 
shown in the plot and is found to be 0.034. 



essary. The RMS scatter here is 0.034, compared to ~0.04 at the 
AAT Due to the number of uncertainties that can be introduced 
in the reduction procedure, such as the scattered light removal, 
blaze extraction, flatfielding etc, we believe that this RMS error 
is the best overall estimate of the errors associated with each in- 
dividual activity measurement. When we compare the final Hall 
et al. log7?[jj. values for the two calibration stars HD1835 and 
HD 10700 against the values presented here, we see good agree- 
ment between the two within the uncertainties expressed i.e. - 
4.41/-4.43 and -4.94/-4.91 respectively. AU final log/^^KFEROS 
indices, along with their associated 5-in dices, J ohnson colours 
and visual magnitudes from Hipparcos dPerrvman et al.,, 19971) 
are shown in Table |4] 

3.2. Metallicities 

When attempting to obtain accurate metallicity abundances 
([M/H]) for solar- type stars some prior knowledge of the stel- 
lar envelope conditions can be useful. Various important stellar 
parameters are input into spectral synthesis codes to model stel- 
lar spectra, these are the effective temperature (Teff), the stellar 
surface gravity (logg), the microturbulence parameter (v^c), the 
macroturbulent velocity (v,nac) and the stellar rotational velocity 
(viot). All of these parameters serve to either alter the line profiles 
and/or strengths of elemental abundances in the stellar envelope 
and can have marked effects on any abundance determination. 

3.3. Teff Determination 

Various tables consisting of numerous stellar parameters and 
calibrations have been published in the literature. For example, 
Stromgren w/jvy-pho tome try for a host of stars can be found in 
Hauck & Mermilliod ( 1998) and using calibrations from Olse^ 
d 19841) the effective temperature s and surface g ravities can be 
found. Also the catalogue of Nordstrom et al.l (12004) has an 
extensive list of effecti ve temperatures for over 16,000 stars. 
iNordstrom et al.l (|2004|) has also published the Stromgren de- 
rived [Fe/H] values for all stars in their catalogue and, along 
with the calibrations o f Haywood (2002) for solar-type stars 
in the iHauck & Mermilliod (il998i) catalogue, a large library of 
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Stromgren [Fe/H] can be compiled. A large number of stars 
in this study have measured Stromgren indices. Such measure- 
ments only allow for an estimate of the overall metallicity abun- 
dance. Spectral analysis is a more accurate technique for the de- 
termination of the overall metallicity abundance, albeit at greater 
labour, as the information contained in each line is measured in- 
dividually and then combined. Whereas using photometry means 
losing all the information on each individual element and relying 
on the stability and accuracy of the photometry and also the cali- 
bration of the scale with stellar parameters. Hence, all stars with 
Stromgren indices were kept in this observing program but pri- 
ority was given to stars with no information at all. Also as most 
stars in this project did not have measured Stromgren indices a 
different method was needed to obtain the stellar temperatures. 

Recently a common method of temperature determination 
is by use of fitting the stellar spectra themselves (e.g., se e 



ii 



ISantos et al.|[2003l: ISantos et al.|[2004t IValenti & FischShoOSh . 
However this method can be quite computer intensive and quan- 
tifying the errors can be quite difficult. A simpler method for ob- 
taining the effective temperature is by using photometric colour- 
effective temperature relations (e.g. . Smith 1995) (equ ation 8.9), 



which was recently employed by iBond et alj (|2006|) to deter 



mine the effective temperatures of st ars on a full target list. The 
Infrared Flux Method (IRFM) (Blac kwell & Lvnas-Gravlll994 
was the calibration chosen to generate the effective temperatures 
using the Hipparcos V magnitude and the 2MASS ^s magnitude 
of all these stars. Since the 2MASS ^s is a shorter band measure- 
ment than the typical Johnson /T-band, a magnitude correction 
taken from Carpenter (2001) was used to correct the K^ magni- 
tude for application to the IRFM. The calibration from Blackwell 
& Lynas-Gray was used to generate the final effective tempera- 
tures for all stars. The V - K colours have been found to be 
more reliable as a temperature indicator due to less metallicity 
independence, reduced line blanke ting, the colour var iation is 
less susceptible to surface gravity dAlonso et al.llI996l) and the 
colours span a wider range of values compared to Z? - V colours, 
giving a tighter sequence. 



3.3.1. Surface Gravity 

The surface gravity (log^) of stellar objects is another parameter 
that can alter stellar line strengths, therefore altering the final 
abundances measured. For each of the stars in this list the x^ for 
each chosen absorption lines was measured for a range of logo's, 
3.5-5.0 dex in steps of 0.5. This allowed a wide parameter space 
to be probed and to test if the final metallicities were drastically 
affected by the surface gravity. It was found from x^ fitting that 
the metallicity of these objects are only weakly dependent on 
the surface gravity and by interpolating between the points in 
;^f^-space the error introduced is minimal. 

However, the analysis in Section 5.3 requires some knowl- 
edge of the stellar surface gravity to use as a proxy for evolu- 
tionary status and this was done by searching the stellar spectra 
for a large region that would produce significant enough gravity 
sensitivity within our course grid of gravity models. The region 
of spectra chosen was centered at 6200A and had a width of 
200A. By extracting this region from all stellar spectra and com- 
paring it to the same region extracted from the models, for each 
of the gravity values and best fit temperatures and metallicities, 
then the x^ minimisation could be used to indicate the evolu- 
tionary status of the star. Since this was across a course grid 
in surface gravity space only an indication could be made, and 
hence construction of a grid with higher resolution would be re- 



quired to determine the actual surface gravity with any degree of 
certainty. However, a check of a number of objects that overlap 
Valenti & Fischer (2005|) and [Gray (2005) have shown that the 
X^ can determine the minimum and a selection can be performed 
to remove any low surface gravity objects. For instance, the star 
HD27442 has been shown to have a \og{g) of 3.78 by Valenti & 
Fischer and 2.93 in Gray et al. and the;^^^ here minimises at 3.5, 
which is the lowest gravity value in this study. This shows the 
technique can be used to select against too evolved and/or too 
young stars. 

3.3.2. Microturbulent Velocity 

Another parameter that has been introduced to smooth over the 
differences obtained between the predicted and observed equiv- 
alent widths of absorption lines as a function of line strength is 
the microturbulence parameter (vmic)- This arbitrary parameter is 
not necessarily required in 3D convective models of stellar atmo- 
spheres ( Asplund et al. 2000) but has become a standard when 
dealing with one-dimensional stellar spectra to alleviate any dis- 
crepancies in equivalent width measurements. Vmic is used to rep- 
resent uncertainties in a number of parameters like excitation 
energy, oscillator strength, effective temperature, surface gravity 
etc. Valenti et al. find that v^ic and the total stellar abundance 
([M/H]) are partially degenerate therefore they fixed the micro- 
turbulence velocity to 0.85kms ' to minimise errors in their final 
abundance values. In this study the microturbulence parameter 
has already been fixed to 2kms ' as the lKuruca (11993 ) model at- 
mospheres were computed for this value across all the required 
metallicity range and also to create a consistent sample. This 
does not introduce any significant uncertainty since our line list 
selection aimed to target weak lines that are formed deeper in the 
stellar interior and hence occupy the linear part of the curve of 
growth. Such lines are insensitive to changes in the microturbu- 
lent velocity. Hence, we believe that the final x^ estimated error 
represents the uncertainty of this parameter. To fully test this we 
computed 10 values using microturbulences of both Ikm/s and 
2km/s for solar-type stars with a range of metallicities and found 
no offset of the final abundances, with only a small RMS scatter 
+0.02 dex. 



3.3.3. IVlacroturbulent Velocity & Stellar Rotation 

In the first approach the macroturbulent velocity (vmac) parame- 
ter has the same affect on atomic line profiles as that of rotation 
or the instrumental profile, however it is described by a different 
profile than both these parameters. As it's name suggest, Vmac 
is the term that considers the smaller scale motions of the stel- 
lar envelope induced by internal motions within the star. High- 
resolution spectra and detailed imaging of the Sun's surface that 
considers the larger scale motions has r evealed velocity field s 
that vary typically at a few kms ' (e.g. ' Rieutord et al.ll200l1) . 
These motions are generated by photospheric granulation and 
acoustic waves permeating throughout the star due to subsurface 
convection. This motion is also present in the high-resolution 
stellar spectra generated from the integrated flux from a stellar 
source and manifests itself by broadening line profiles. However, 
as the average resolution across the spectrum is ~0.15A, (a ve- 
locity of ~6kms"'), we believe the broadening to be dominated 
by the stellar rotation as velocity fields of a few kms"' are com- 
parable only with the slowest rotation rates. 

As mentioned above, stellar rotation (vsin/) has the effect 
of broadening line profiles without altering the line equivalent 
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width, therefore knowledge of the rotational velocity allows a 
better representation of the line profile and will further reduce 
the errors on the final metallicities. In order to determine the 
vsin/ for each star we selected an iron line that was shown to ex- 
hibit little change in line strength through the range of metallici- 
ties. This line, which is centered at a wavelength of 7445.758A, 
is highlighted by an asterisk in Table [1] and is the strongest line 
in the list. The observed spectra around this line were then com- 
pared with the model spectra, broadened in steps of Ikms ' be- 
tween l-15kms"'. We chose a width of ±0.5A centered on the 
line to compare each line with the models as this also included 
valuable wing information, x^ values were generated at each ve- 
locity step and the minimum represents the final vsin/. The reso- 
lution of the instrument at 7000A is ~6kms" ' , however we found 
that we can discern values of vsin/> 2.5kms"' . 



3.3.4. Model Atmospheres 

To enable highly accurate chemical abundances to be deter- 
mined from this dataset a number of utilities were used to 
generate relevant synthetic spectra to compare with the ob- 
servations. The prog ram chosen t o creat e the s ynthetic spec- 
tra was WITA6 ( Pavlenkoll2000h . iMoord (11956) was used for 
the atomic line list an d damping cons tants, with the ATLAS 9 
mod el atmospheres of Kurucz ( 19931) and oscillator strengths 
from lGurtovenko & KostikI d 19981)" 

The grid of Kurucz models range from 3500K to 8000K with 
surface gravities and metallicities from 3.5 to 5.0 and -1.5 to 
1 .0 respectively, ensuring coverage of almost all possible values 
for these parameters for all objects. A few objects require lower 
metallicity models to accurately measure their values. A param- 
eter that is fixed throughout is the microturbulent velocity. This 
is held fixed at 2kms~' for all the model grids, which will intro- 
duce some error into the accuracy of the final synthetic spectrum 
but after investigation it was found to be small compared with 
other uncertainties. 

Since most of our objects are G and K type stars these 
models are adequate f or metallicity computation, however 
ISmallev & Kupkal ( Il997h have shown that there are errors inher- 
ent in them due to an inadequate convective overshooting algo- 
rithm. Although the errors are small for G and K sta rs, they can 
be significant for hotter F type stars (* Grav et al1l2003 ). To inves- 
tigate this we compute the metallicities for all objects with effec- 
tive temperatures greater than 6000K, with and without convec- 
tive overshooting. The results from the convective overshooting 
On-Off gave a mean of -0.01+0.06 dex for a total of 19 data 
points, which is consistent with an additional scatter of 0.06 dex 
for the hottest stars in this sample. 
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3.3.5. Line Selection 

One of the most critical parts of abundance determination is the 
selection of lines to use in the observed spectrum to obtain a 
clean line profile. A vast number of the lines in a star such as the 
Sun are blended, especially towards the bluer end of the visual 
band (see Fig.|4]upper panel). The lines used in this work have 
been selected because they appear unblended at this spectral res- 
olution (Fig. ID lower panel) and they are also weak in the solar 
spectrum. Weak lines lie on the linear part of the curve of growth, 
therefore they are more sensitive to large changes in the model 
parameters and are better indicators of line profile characteris- 
tics giving more accurate abundances. A line's position on the 
curve of growth will significantly aff'ect all parameters associ- 
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Fig. 4. A spectrum of the star HD38459 towards the bluer end of 
the stars visual band is shown in the upper panel. All of the lines 
in this wavelength region (4074-4078A) are blended. The Fel 
line in the middle of this region (4076. 22A) is clearly blended 
at both sides. The lower panel shows the iron line at 7586. 027A. 
No strong line blends can be discerned from the noise at this 
spectral resolution and if any do exist they are rather marginal. 



ated with the line as the curve of growth determines how the line 
evolves when more absorbers are introduced. For example, opti- 
cally thick lines will saturate more quickly and therefore are not 
good estimators of the line strength and hence abundance status 
of the star. Strong lines also add extra uncertainty through satura- 
tion, microturbule nce and unc ertainties associated with damping 
constants (e.g. seel Gravll2005l) . 

The FEROS spectra used here have a S/N and resolution 
sufficiently high enough to allow a large enough number of ap- 
parently unblended lines to be used in this analysis. The lower 
wavelength cutoff chosen was at a /I > 5400A this was due to 
the increased blending factor mentioned above (Fig.|4|, however 
any lines below 6500A were very carefully selected to ensure 
they were isolated from other possible blends. The wide cover- 
age region obtained using FEROS (-3500 - 9200A' Kaufer et al.1 
1999) allows one to select a large enough number of unblended 
iron lines. Solar spectra from FEROS, taken as part of the stan- 
dard ESO Calibration plan, was used as a template to select lines 
that would appear in all other spectra. Also, no telluric lines were 
seen to blend with these lines. All lines used for abundance de- 
termination are shown in Table. [1] 

3.3.6. Oscillator Strengths 

One of the larger sources of error comes from uncertainties 
in the oscillator strength values. The oscillator strengths cho - 
sen were extracted from the list of iGurtovenko & KostikI (Il998l) . 
Gurtovenko & Kostik measure the oscillator strength (gf) using 
measurements of the equivalent widths, central intensities along 
with the synthesis of the whole Fraunhofer spectrum of the spe- 
cific chemical element. They define two separate measurements 
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Table 1. Fe lines used in the abundance measurements are tabulated with their wavelengths, oscillator strengths and excitation 
energies. The iron line highlighted by the asterisk was used in our determination of v sin /. 



Element 


^(A) 


gf 


E (keV) 


Fei 


5806.723 


1.096E-01 


4.610 


Fei 


5852.215 


5.623E-02 


4.550 


Fei 


5855.082 


2.239E-02 


4.610 


Fei 


5856.086 


2.291E-02 


4.290 


Fei 


6027.051 


5.888E-02 


4.070 


Fei 


6151.621 


4.266E-04 


2.180 


Fei 


6159.379 


1.148E-02 


4.610 


Fei 


6165.359 


2.754E-02 


4.140 


Fei 


6173.340 


1.259E-03 


2.220 


Fei 


6229.227 


9.120E-04 


2.840 


Fei 


6608.027 


9.120E-05 


2.280 


Fei 


6627.547 


2.570E-02 


4.550 


Fei 


6703.566 


8.318E-04 


2.760 


Fei 


6725.359 


5.370E-03 


4.100 


Fei 


6750.152 


2.344E-03 


2.420 


Fei 


6810.266 


8.710E-02 


4.610 


Fei 


7158.477 


1.148E-03 


3.650 


Fei 


7189.152 


1.585E-03 


3.070 


Fez* 


7445.758 


6.761E-01 


4.260 


Fei 


7491.656 


8.710E-02 


4.300 


Fei 


7568.906 


1.288E-01 


4.280 


Fei 


7586.023 


7.079E-01 


4.310 


Fei 


7751.109 


1.514E-01 


4.990 


Fei 


7780.562 


1.148E+00 


4.470 


Fei 


7807.914 


2.754E-01 


4.990 


Fei 


7832.207 


1.259E+00 


4.430 


Fei 


8047.625 


2.239E-05 


0.860 


Fei 


8239.137 


3.715E-04 


2.420 



of gf using separately the equivalent widths of the lines and the 
central line depths to allow them to estimate the accuracy of the 
values and to generate estimates of the uncertainties in the data. 
They found that the internal accuracy of the final lo^igf) values 
to typically be ±0.07 dex. Recendy iBigot & TheveninI (|2006|) 
have determined the oscillator strengths for Fei and Sil in the 
Gaia spectral window (~8480-8750A) and found accuracies of 
typically < 0.1 dex with laboratory results, whereas the oscil- 
lator strength accuracies found in databases such as the Vienna 
Atomic Line Database can be much lower. 



3.3.8. Continuum Normalisation 

No flux calibration was performed which left behind a gradient 
in the spectral profile across the measured wavelength region. 
To attempt to detect the true continuum level around each line 
a cubic-spline fit procedure was performed using the program 
CONTINUUM part of the iraf library of routines. A global fit was 
used across the entire wavelength region. By employing tight 
constraints on the upper and lower limit rejection thresholds for 
the fitted data points a good continuum measurement was made 
for the entire spectra. 



3.3.7. WITA Program 



The latest version in the WITA series of programs was used to 
generate the grid of synthetic spectra. Th e WITA series w ere de- 
veloped as part of the ABEL6 complex ('P avlenkol ll99l'). which 
was composed partly from Kurucz ( 1979) ATLAS subprograms. 
The WITA6 program (Pavlenko 200 0i) itself is a m odified ver- 
sion of the WITA2 program (Pa vlenko et al.|[l995l) . All compu- 
tations performed by WITA6 are employed in the classical way, 
using Local Thermodynamical Equilibrium (LTE) and assum- 
ing hydrodynamic equilibrium for one-dimensional model atmo- 
spheres. The code allows for the extrapolation of the model at- 
mosphere. The ionisation-dissociation equilibria equations sys- 
tem has been solved for various temperature structures assuming 
LTE. A Voigt profile was adopted for single absorption lines, tak- 
ing into consideration all line broadening sources, such as nat- 
ural broadening, thermal broadening, pressure broadening and 
resonance broadening. 



3.4. x^ Fitting 

Once all the parameters for each object have been fit and iron 
lines selected the final step is to generate the desired [Fe/H] 
abundance. An IDL code has been developed that will employ 
several different steps to compare the lines in the real spectrum 
with the corresponding lines in the synthetic spectrum and gener- 
ate ax^ value for each line in each star These ;i'^'s will minimise 
when the synthetic spectrum best replicates the real spectrum, 
allowing [Fe/H] abundances to be determined. 

The first step in the program takes the generated tempera- 
ture values for each star and then selects the spectra that have 
the nearest temperatures above and below the generated temper- 
ature in the grid. The Kurucz grid has temperature steps of 250K, 
therefore to better replicate the real spectra the upper and lower 
spectra are interpolated to an improved temperature scale. It is 
at this stage when the synthetic spectrum is broadened for stel- 
lar rotation and the instrumental profile. The rotation rates are 
generated following the procedure described above. The instru- 
mental profile is modeled by convolution with a Gaussian profile 
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Fig. 5. The x^ fit for HD1835 in metallicity space. The points 
represent the means of four ;^'^ values for different surface grav- 
ities at each of these metallicities. The solid line is the second- 
order polynomial fit to the points, with the final metallicity rep- 
resenting the minimum of this function. 

of width determined by the resolution of FEROS. The FEROS 
resolution is ~48'000, which at these wavelengths corresponds 
to a Gaussian of width ~0.15A however since the instrumental 
profile changes with wavelength (e.g. telluric lines in the far red 
of the chip have widths >0.2A) we also use the models to better 
define an optimum broadening width. 



X- — --, — 



(2) 



After the model spectra have been prepared for comparison 
with the real spectrum the program selects each individual line 
that has been selected for comparison and employs a x^ fit be- 
tween the real and model spectra. Equation. |2] shows the equa- 
tion employed for this fitting procedure, where O is the observed 
spectra, M is the modeled spectra and o" is the total errors in 
both. This x^ equation is applied to each line by selecting the 
line center and comparing a region +0.2A from the center The 
0.4A window was selected by eye after scrutinising all lines and 
finding the optimum window for these lines in the FEROS spec- 
tra. As this window is only an optimal window and will include 
more continuum in some lines than others, the median of all the 
X^ values is used. This helps remove any outlying values due 
to bad lines, any possible small edge blending and extraction 
errors. The median is generated for each metallicity value four 
times, representing the four different log(g)'s for each metallic- 
ity. Figure |5] shows the fitting output for the star HD1835. The 
crosses are the medians of the four separate log(g) fits and the 
solid line represents the bestfit second-order polynomial to the 
points. The minimum of this function represents the stars [Fe/H] 
abundance, which for th is star is 0.21+0.03, agreeing with the 
IValenti & Fischerf ( l2005h value of 0.25+0.03. This entire proce- 
dure is applied to all targets. On a standard Pentium 4 processor 
with 1 Gigabyte RAM, it takes around 2 days to complete. 

3AA. x^ Uncertainties 

To determine the errors in the final abundance measurements the 
fitted ;(f^ parabola was used to interpolate the 1,2 and 3cr values 
using a table of;t'^ values. Figure |6] shows the;t'^ space for all 
three of these error ranges. The solid lines bound the Icr error 




0.2 0.3 0.4 

[Fe/H] 

Fig. 6. The;^'^ fit across all metallicity values for the star HD1835 
with the Icr (solid bounds), 2cr (dashed bounds) and 3cr (dot- 
dashed bounds) error estimates shown. These errors ranges were 
selected by interpolating a table of ;\f^ steps for a two parame- 
ter fit i.e. surface gravity and metallicity. The central solid line 
marks the minimum of the fitted function and represents the final 
abundance value derived for this star. 



region, with the central solid line marking the minimum of the 
function or the final abundance measurement. The dashed and 
dot-dashed regions mark the 2 and 3cr parameter space. Since 
the function is fit to both the \og(g) and the iron abundance a 
two parameter ;t'^ error step was performed. The Icr errors have 
a step of 1 . 155 from the minimum value, with the 2 and 3cr steps 
equaling 3.118 and 6.396 respectively. This error analysis is only 
valid assuming one has a handle on the uncertainty, however as 
shown before their are a number of uncertainty sources that are 
difficult to quantify, (i.e. atomic oscillator strengths), therefore 
the final x^ curve must represent this uncertainty. To accomplish 
this the reduced Chi-squared (x^) was also determined by divid- 
ing the x^ by the number of degrees of freedom. The degrees of 
freedom were found by subtracting the number of fitted param- 
eters from the number of data points for each iron line. Ideally 
the xj should equal, or be close to 1 for a good fit with all uncer- 
tainty sources well modeled. However, since many uncertainties 
are poorly constrained, the x^ function was normalised by the 
difference between the minimum of the x^ function and the op- 
timum (i.e. 1). This has the same effect as randomly adding per- 
centage errors into the original x^ equation (Eq. |2]) to broaden 
the fitting function and increase the final uncertainty estimates 
for all abundances. 



4. Comparison With Previous Work 

4.1. Activities 

To make judgments and comparisons of this data to other pub- 
lished work we have plotted the individual activity values against 
those of three other catalogues. Fi gure 111 shows the compar- 
isons o f these values to those of iHenrv et al. (1996) (upper 
panel), IWrightl (l2004al) (middle panel) and iGrav et alj (|2006|) 
(lower panel). The comparison to that of the Henry et al. val- 
ues shows a scatter around the 1 : 1 relation (solid line) with no 
apparent offsets. Considering they also used calibrators from 
Duncan et al. (1991) one would expect there to be a good re- 
lationship. However, the Gray et al. work makes use of low- 
resolution spectra and also uses a different setup from that em- 
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Fig. 7. A direct cornparis on of logR^p- values from this work 
againstlHenrv et al.l(ll996h (upper).IWnght (2004a) (middle) and 
iGrav et al] (l2006h (lower). The plots show good agreement with 
Henry et al. and Wright et al. and both have RMS scatters of 
~0.08 in log space. However, there appears a small offset with 
the results of Gray et al., with the values here slightly higher by 
around 0. 1 dex. The solid line through the data represent a 1 : 1 
relationship. 



ployed here. They centered 4A square bandpasses on the Can 
HK line cores instead of triangular bandpasses of width 1 .09A 
which is used at Mt. Wilson and also used here. This will con- 
tribute to an increase in the scatter between the data and also to 
non-linear calibrations onto the MW system. The RMS scatter in 
this data is found to be ~0.08. A Student-T test was performed 
on all three of these distributions to test if they have the same true 
means. The test revealed an 83% probability that the data from 
here and that of Henry et al. are drawn from the same mean, 
which given the intrinsic variability associated with cool dwarf 
stars and the limited number of data points is noteworthy. 

The middle panel in Fig. [T] shows the data from Wright 
( l2004ab against the values found here. There appears a good fit to 
this data around the 1:1 relation with only two outliers. This is a 
nice fit when considering the intrinsic scatter inherent in the val- 
ues themselves, which can alter by 0.15 in log/?[jj^ over a period 
of a few years. The RMS scatter here is ~0.08, the same as the 
scatter of Henry et al. However, the two outlying stars (HD8038 
and HD 150437) cause this unweighted RMS to deviate from the 
true scatter. Once these two stars have been removed the RMS 
scatter is found to be 0.06, highlighting the close relationship 
that is found using high-resolution spectra and a setup similar 
to that of the MW project. This is highlighted by the fact that 
the T test for this data reveal that they are both drawn from the 
same mean at a probability of 88%. Again this is high given the 
limited amount of data points to compare against. 

A recent large survey in th e southern hemisphere has been 
conducted bv lGrav et al.l (l200 6l) connected to the NStars Project. 
The lower panel in Fig. [T] compares our measurements to this 
data. There is some scatter around the 1 : 1 relation, however the 
majority of objects lie below this fit indicating a possible sys- 
tematic offset in the data. Neglecting the stars above the 1 : 1 fit 
there is an offset of around -0. 1 in XogR'^^^ that is not present in 
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Fig. 8. A compar ison of the \op i '^ y activity indices from 
iHenrv et al.1 d 19961) and lGrav et"aLr(l2006l) . There is an offset in 
the data, with the values of Gray et al. around 0. 1 dex lower than 
those of Henry et al. 



both other studies. A straight line fit to describe the offset has 
a gradient of 0.938+0.006 and off^set of -0.389+0.362, with an 
RMS scatter of 0.096. This can be used to transform these re- 
sults to the scale of Gray et al. The significance of the offset is 
at the 2cr level as the T test for this comparison gave a proba- 
bility of 96% that both distributions were drawn from different 
true means. We also performed the log/?J^j^ conversion includ- 
ing all the stars in common with Gray et al., however once this 
calibration has been performed, both Henry et al. and Wright 
are offset from this work by around +0.1 dex. The point can be 
highlighted by plotting Gray et al. against the work of Henry et 
al. as both studies used similar resolutions and setups (Fig. [HJ. 
It can be seen from this comparison that the Gray et al. dataset 
is around 0. 1 dex lower in log/?[jj^, which is difficult to explain 
by intrinsic scatter. Indeed, the T test revealed a probability of 
98% that both have been drawn from different true means. The 
offsets could arise due to the different calibrators used by this 
and the Henry et al. study against those of Gray et al. Gray et al. 
used calibrating stars from Bali unas et al.l (Il995h compared to 
the Duncan et al. values, however as no offset appears between 
this study and Wright, (who also calibrated from Baliunas et al.), 
this seems unlikely. It may simply arise from a combination of 
the calibrators and the wide bandpasses used to extract the ac- 
tivities as the width, shape and resolution in the line cores can 
serve to mask any emission feature and hence mask the true ac- 
tivity of a star. Note also that Henry et al. performed non-linear 
calibrations onto the MW system, whereas both this work and 
Gray et al. use linear calibrations and this may explain some of 
the offsets seen, especially towards the more active end of the 
calibration. Whatever the cause it is clear a more thorough look 
at all these studies together, along with long term variability data 
is needed to extract a good list of intrinsically stable stars to use 
as calibrators for similar studies in the southern hemisphere. 

A.2. Metallicity 

An integral part of any abundance determination of a large sam- 
ple of solar-type stars is to compare the final values with those 
published in other work. We have chosen three of the most re- 
cent large bodies of published metallicity values in the south- 
ern hemisphere to compare against our determined abundances. 
The th r ee works we have compa red against are Valenti & Fischej 
( l2005l) . lNordstrom et alj(l2004 and lGrav et al.l ( l2006l) . The most 
appropriate work to compare the results found here with are 
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Fig. 9. A direct comparison of effective temperatures (left panels) and metallicity values (r ight panels) deterinined h ere and those 
published in recent work. The upper plots are a comparison with the values published in IValenti & FischeJ (l2005h . which were 
derived using spectral line fitting similar to the analysis in this project for metallicities but differing in effective temperature de- 
termination. The middle plots are a comparison with the values of Nordstrom et al. (2004), which were derived from Stromgren 
photometry. The lower plots compare against the values published in lGrav et al.l ((2006), derived by fitting to low-resolution spectra. 
The effective temperatures agree well with those of Valenti & Fischer, with only two outliers that are variable stars. There appears 
to be offsets between this work and both Nordstrom et al. and Gray et al. The metallicity plots agree well with all three of the 
published work, with the lowest RMS scatter found against the work of Valenti & Fischer due to the similarity in the methodology. 
However, a scatter of only ~0.1 was also found compared to the other methods. The solid lines through all six plots represent the 
1:1 relationship. 



those of Valenti & Fischer as the methodologies are similar 
They determine their final abundances by directly comparing ob- 
served high-resolution, high-S/N echelle spectra against a grid of 
synthetic model spectra generated using Kurucz (1992) model 
atmospheres. This comparison is shown for both the effective 
temperatures and the metallicities in the upper panels (left and 
right) of Figure |9] The TEFF.Vaienti and [Fe/H]vaienti are the ef- 
fective temperatures and iron abundances by Valenti & Fischer 
and the Teff.feros and [Fc/HIferos are similar but for the values 
found in this project. The solid lines represent a 1 : 1 relationship 
and it is clear that the values in both panels are in good agree- 
ment. The effective temperatures (left panel) in Valenti & Fischer 
were derived by a spectral fitting procedure, which is different to 
our photometrically derived IRFM values, however the data are 
in good agreement with only two outliers that are significantly 
below the 1 : 1 relationship. Both outliers have been found to ex- 
hibit some level of variability by Hipparcos and hence the pho- 
tometric method can produce spurious results. The RMS scat- 
ter around a linear fit to this data is +132K, which reduces to 
+80K once the two variable outliers have been removed. The 
Student-T test reveals a 85% probability that both these have 
similar means, increasing to 91% when the outliers have been 
removed. The metallicity comparison (right panel) shows a tight 
relationship between that data. This is confirmed by the scat- 
ter of only +0.09, reducing to ±0.05 once the two outliers have 
been removed. The T-test revealed a 96% probability that both 
distributions have the same means, reducing to 95% when the 
outliers have been removed. However, if we consider all objects 
with metallicities >0. Idex from both Valenti & Fischer and here, 
then the Valenti & Fischer values are found to be more metal-rich 
by -H0.025dex. 

The middle and lower panels in Figure |9] show compar- 
isons against the two other southern samples, however the 



methodologies were different to that employed here. The mid- 
dle panel shows the co mparison with the values published in 
Nordstro m et al.1 (|2004^ ([M/H]Noi-dstrom)- Nordstrom et al. de- 
termined the metallicity of ~ 14000 stars using relationships be- 
tween the overall stellar metallicity and their Stromgren colours. 
Stromgren photometry has proven a useful technique in estimat- 
ing a star's metallicity, however the calibrations tend to have 
larger systematic uncertainty than spectroscopically determined 
abundances due to a lack of calibrators, particul arly towards 
the reddest stars studied (e.g. iTwarog et al.l l2002h . Nordstrom 
et al. created their own calibrations following techniques previ- 
ously employed by othe r authors (e.g., Schuster & Nissen 1989f 
lEdvardsson et al.ll 19931) and they find a scatter of 0. 12 for G and 
K-type stars, reducing to 0.1 for F-types. It must be noted that 
the relations determined are compared to the iron abundances 
([Fe/H]) of the calibrator stars since iron is an accurate tracer 
of the overall stellar metallicity abundance. From Fig. |9] (mid- 
dle left) it can be seen that the scatter in effective temperature 
between this work and Nordstrom et al. is lower than both the 
other two works. The measured RMS scatter is 73K, however a 
step change of +50K seems apparent at an effective temperature 
of 5250K. For the whole sample the T-test reveals a 82% proba- 
bility that both these distributions have similar means. However, 
for the lower and upper temperature cut the probabilities are 5% 
and 28% respectively. These offsets mainly arise due to the dif- 
ferent methodologies used to generate the effective temperatures 
i.e. Stromgren versus Johnson filters. The metallicity compari- 
son (middle right) has a scatter of 0. 1 3 dex and a T-test probabil- 
ity of 100% that both are drawn from the same mean. This level 
of agreement between the two data show that with a sufficiently 
large comparison dataset, the values in this work are as accurate 
as any measurements currently in the literature. 
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As mentioned in Section 4.1, the Gray et al. work is part 
of the NStars Project and their methodology employed fitting 
models directly to a large spectral region of low -resolution spec- 
tra, whilst also fi tting to the observed and synthetic fluxes (see 
Gray et al. 20031) . The synthetic spectra were generated using 
Kurucj (1 19931) model atmospheres and the spectral synthesis 
program spectrum . A Chi-squared minimisation technique was 
used, in conjunction with a simplex interpolation to generate a 
global fit to the observed spectra. Again Gray et al. provides a 
useful test as their method made use of stellar spectra to fit mod- 
els directly. They have also compared t heir overall [M/HI values 
to the [Fe/H] abundances measured bv lValenti & Fischer (l2005h 
and found a scatter of only 0.09 dex, however they found an off- 
set of -0.07 dex. From the lower right plot in Fig. |9]it can again 
be seen that the derived metallicity values are in good agree- 
ment with those of Gray et al. ([M/H]Gray)- The scatter is found 
to be 0.09 dex, which compares well with the scatter found be- 
tween Gray et al. and Valenti & Fischer. There may be an off'set 
of -0.01 dex between these samples, however since this is be- 
low the general scatter in the data, this is deemed negligible. 
This is confirmed by the T-test which reveals a probability of 
81%, indicating that the offset is not significant. The lower left 
plot in Fig. |9] shows the comparison with the Gray et al. effec- 
tive temperatures and a good coiTelation is found between the 
two datasets. The scatter is found to be +90K and no signifi- 
cant offset was found, giving rise to a probability of 93% that 
both datasets have the same means. Since Gray et al. had shown 
their values to be in agreement with those of the IRFM, a good 
correlation was to be expected. By comparing FEROS values to 
three independent studies that employed three different method- 
ologies it is clear that these values provide a robust assessment 
of the effective temperatures and iron abundances of all stars in 
this sample. 



5. Results & Discussion 

5.1. Activity Analysis 

The distribution of chromospheric activity for solar-type stars 
has been shown to follow a b imodal distribution with both an 
active and inactive pe ak (e.g. iDuncan et al.l 1 1 99 It iHenrv et al] 
Il996l:l&av et al.l2006t J06). The upper panel in FigurefTOlshows 
the distribution of stellar logR'^^ in this study and confirms the 
bimodal peaks found before. The plot has been binned up with 
bin widths of 0.05 in logR'^^ and both peaks are easily identifi- 
able. A bimodal gaussian has been fit to the data and is modelled 
by the dashed line in the plot. The peaks of this distribution are 
centered at -4.52 (active) and -5.00 (inactive) and have standard 
deviations of 0.13 and 0.10 respectively. There also are a few 
stars out in the tails of the distribution. HD27442 and HD 147873 
are located out in the inactive wing and have logR'^^^ indices of 
-5.39 and -5.42 respectively. 

The lower plot in FigurefTOlshows the residuals of the dou- 
ble Gaussian fit to the distribution. The stars located out in the 
active tail of the fit are single objects and therefore the residual 
errors bars will place them in agreement with the fit as the un- 
certainty is approximated using Poisson root-N, statistics. From 
these statistics alone it is hard to test whether or not these stars 
are undergoing an active or inactive phase of their evolution. 
Almost all the bins agree well with the fit within the uncertain- 
ties shown except for one which is over 2cr away from the fit. It 
will be shown later in the chapter that this bin, which is centered 
at -5.175, contains mostly evolv ed stars. It is inter esting to note 
that this result was also found bv lGrav et al.l (l2003h at around the 
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Fig. 10. The upper plot shows the logR'^^ distribution of all 
bright stars in the southern hemisphere in this study and the in- 
active and active peaks are confirmed. The overall distribution 
is i n agreement with t he bimodal fits from iHenrv et alj (1 19961) 
and Grav et aTl ( l2006l) in the southern hemisphere, however a 
third component may exist. The best-fit to the data employs a 
bimodal distribution and this is overplotted by the dashed line. 
The means of the Gaussians are -4.52 (active) and -5.00 (inac- 
tive) with standard deviations of 0. 13 and 0. 10 respectively. The 
lower plot shows the residuals of the fit along with their asso- 
ciated error bars. Most of the bins are in good agreement with 
the fit, with only one bin over 2cr away from the fit. This bin 
contains mostly stars evolving off the main sequence. 



same inactive bin of their distribution. They concluded that this 
may partly be caused by a larg e number of Maunder Minimum 
stars. However, Wrighr ?2004tf) found that most stars announced 
as in Maunder Minimum phases were in fact stars evolving off 
the main sequence and hence the activity methodology used to 
measure their indices are not entirely adequate for comparison 
with main sequence stars. A trimodal fit was also applied to the 
data to test if the fit was significantly improved. The x^ for the 
bimodal distribution was 30 with 23 degrees of freedom, com- 
pared with the value of 21 with 20 degrees of freedom for the 
trimodal distribution. Therefore the reduced x^ values are 1 .29 
and 1.06 respectively. Thus, it is not clear that a significant im- 
provement is made by assuming a third peak exists in the data. 

5.2. Metaiiicity Analysis 

To assess the validity of the metallicity abundances for each 
of the target stars a check of line excitation energy and wave- 
length dependency is performed. The final abundance for each 
line is computed following the same procedure as above. This 
abundance is then subtracted from each star's total abundance 
and plotted against both excitation energy and wavelength for 
the given line. Fig. [TT| shows both the excitation energy (up- 
per panel) and wavelength (lower panel) abundances for each of 
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Fig. 11. A check of the dependency of the final abundance 
measurements on both absorption line excitation energy (up- 
per panel) and wavelength (lower panel) for the star HD28013. 
These plots are typical of the whole sample. This shows that the 
final abundance values are not significantly dependent on both 
excitation energy or wavelength. The dashed lines represent best 
fit Unear trends to the data. The solid line marks the zero points. 



the lines chosen here in the solar analogue (Spectral Type G5V; 
[Fe/H]=0.01) HD28013. The majority of data points are spread 
within +0.1 dex from the final abundance value in both plots, 
which is represented by the solid line, and it is these clustered 
points that give rise to the measured abundance value. A linear 
least-squares fit was applied to the data and these are represented 
by the dashed lines in both plots. It is clear that the gradient 
of these lines are very small, indicating that there is very little 
dependence of these abundance measurements on excitation en- 
ergy or wavelength. These plots are typical for the whole sample 
of objects and indicate that the computed [Fe/H] abundances are 
reliable and not a function of the selected lines. 

The overall distribution of metallicities for all stars in this 
sample is shown in the upper panel of Figure[T2l The peak of the 
distribution resides at ~0. 0-0.1 dex, however it is hard to draw 
any inferences about the distribution of metallicities in the solar 
neighbourhood as there are specific biases in the data. For in- 
sta nce, a subsection of the data had Stromgren colours measured 
by iHauck & MermilliodI (Il998h and using calibrations found in 
iHavwoodl (I2OO2I) we were able to obtain a first approximation 
of the overall metallicity for these objects. We then selected a 
subset of the most metal-rich of these objects to follow-up with 
spectroscopy to confirm their metal-rich status and further boost 
the numbers of metal-rich objects that will constitute the final 



target Ust. However, as explained in Section 3.3, the primary list 
of targets were objects with no known high-resolution metallic- 
ity measurements, to further enhance the database of fundamen- 
tal stellar parameters in the southern hemisphere. Also shown 
in this figure is the total number of e xpected stars following th e 
probability of hosting a planet from iFischer & Valentil (|2005|) . 
As mentioned above, Fischer & Valenti have found that the prob- 
ability of hosting a planet follows a power law increasing to 
25% for stars with metallicities over -1-0.3 dex. By assuming this 
power law for the sample we have estimated the total number 
of planets expected to be found in this metal-rich sample to be 
~12 and is shown in the upper right of the figure. Note this is for 
all objects not currently on any large planet search programs. Of 
the 137 currently discovered planets around metal-rich stars, 15 
have been found to transit their host star, giving rise to a frac- 
tion of ~11% that any detected planet around a metal -rich star 
will transit. Therefore, we expect there to be at least one tran- 
siting planet in our metal-rich sample that can be followed-up 
with high resolution spectroscopy to probe the atmospheres of 
exoplanets. It must be noted that this only estimates the num- 
ber of hot Jupiter-type planets down to a mass of 0.5Mj and a 
radial-velocity precision of 40ms"'. Any additional hot Saturns 
or Neptunes will not be covered by this estimate and therefore 
the true number of detectable planets should be larger than this. 

From the upper part of Figure [12] it is shown that the distri- 
bution of stars in this biased sample mostly followed a Gaussian 
distribution with an extended metal-poor tail. The other panels in 
Fig.[T2]show three metallicity distributions from the surveys we 
have compared to in Section 4.2. The upper middle panel show s 
the distribution of stars taken from Valenti & Fisched (l2005h . 
This sample contains stars mostly from the northern hemisphere, 
how ever as the northern and southern distributions are similar 
(e.g. iGrav et"ani2003t [Nordstrom et al.ll2004l: iGrav e"tani2006h 
and the methodology used was similar to the methodology used 
here, it can provide a useful comparison sample. The peak of 
this data agrees well with this work (upper panel) as it lies be- 
tween 0-0.1 dex, slightly metal-rich. This is largely again due 
to a bias in their sample selection where they selected all stars 
on current large scale planet search surveys, which itself intro- 
duces a number of metal-rich stars as most surveys now have 
metal-rich programs. Also the extended metal-poor tail is appar- 
ent in this data. Howeve r, when co mpared with the lo wer middle 
(Nordstr om et al.ll2004t) and lower dGrav et alj|2006l) panels the 
distributions differ in where the peak lies. Both these samples are 
taken solely from the southern hemisphere (dec <0°) and both 
have peaks below solar metallicity, with the largest sample from 
Nordstom et al. peaking between -0.1 and -0.2. These samples 
are both unbiased samples as their selection criteria is volume 
limited out to 40pc. Nordstrom et al. selected a magnitude lim- 
ited sample of stars with Stromgren photometry that is complete 
for all F and G stars out to 40pc and is kinematically unbiased, 
whereas the Gray et al sample is as yet incomplete but also fo- 
cused on all stars earlier than MO in the Hipparcos catalogue. 
Both these studies should better represent the overall distribu- 
tion of stars in the galactic disk. It is clear that the Nordstrom et 
al. distribution follows a Gaussian profile with only a very small 
indication of a metal-poor tail, yet the Gray et al. profile shows 
a clear extended tail. The extended tail in the distributions arise 
mostly from stars from both the t hick disk and the galactic halo 
and possibly from active K stars dGrav et al.ll2006l) . 

Figure [T3] shows the distribution of all uncertainties for the 
metallicity fits to these stars in this sample. The uncertainties 
appear to follow a normal distribution peaking at +0.065 dex, 
which agrees well with the scatter found in the metallicities 
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Fig. 12. Four plots showing the distributions of metallicities 
taken fr om this work (upper), V alenti & Fischer ( 20 05) (upper 
rniddle ) jNordstrom et alJ (l2004l) (lower middle) and iGray et al.l 
(I2OO6I) (lower). Since all four distributions have significantly dif- 
ferent sample sizes all were normalised to the peak of the distri- 
bution, with the sample sizes indicated on the plots. The labels 
also represent the constituent samples as the upper, lower middle 
and lower panels represent the southern hemisphere, whereas the 
upper middle panel is mostly the northern hemisphere but was 
included as this dataset represents the comparison sample with 
lowest internal errors. The upper and upper middle distributions 
peak around 0-0. 1 dex due to the metal-rich bias inherent in the 
sample selection, whereas both the middle and bottom panels 
should provide a better representation of the overall metallicity 
distribution in the southern hemisphere. 



when compared to other works in the literature. However, an ex- 
tended tail is found at larger uncertainty values. These are mostly 
comprised of low S/N stars that were observed through cloud 
cover. The low S/N means there are less photons making up each 
Fe line and the line profiles are not well behaved and are harder 
to model. 



5.3. Evolutionary Phase 

To determine if any of the stars in this project are in a Maunder 
Minimum phase of their evolut ion we have tried to ascer tain 
their current evolutionary status. iBaliunas & Jastrowl d 19901) de- 
termined that 30% of the MW stars were in a Maunder min- 
imum phase, this was based purely on their low activity sta- 
tu s. This estimate was rev ised to a lower percentage of 10-15% 
by ISaar & BaliunasI d 19921). which has since been cornplimented 
by oth er studies (iHenrv et al.ll 19961; ISaar et al.ll 19981 iGrav et al.1 
l2003h . However, Wright! ( l2004b l) has shown that these stars 
are all either evolved, so the application of the activity index 
does not apply, or are metal-poor, which alters the CaHK line 
wings for such stars so their indices are unreliable. Note that 
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Fig. 13. The distribution of l-cr uncertainties from the chi- 
squared fitting for all stars in this project. The distribution peaks 
at +0.065 dex, with the majority of values spread between 0.03- 
0.10. There are also a number of stars with larger uncertainties, 
which were mostly the stars that were observed through cloud 
cover and have lower S/N. 

activity indices are calibrated for stars thought to be on the 
main sequence and hence the further from the main sequence 
a star evolves the more unreliable the index becomes. Indeed it 
may be necessary to re-determine the activity indices for more 
evolved stars. What is clear is that the Sun went through a sim- 
ilar activity minimum in the 17* Century and it has been es- 
timated that the solar ^mw index decreas ed to a mean value 
of around 0. 145 (IBaliunas & Jastrowll 19901) (currently the mean 
is ^0=0. 165). By applying this determination to the analysis 
of brightness changes due to activity, Zhang et al. ( 1994) have 
determined that the Sun was 0.2%-0.6% dimmer during the 
Maunder Minimum. Stars in this study may be in a similar evo- 
lutionary phase. Figure [14] shows the distribution of stars in this 
survey on a colour-magnitude diagram. The B - V colours and 
absolu te magnitudes were e xtracted from the Hipparcos cata- 
logue (iPerrvman et alJlT997b . The main sequence is clear in this 
plot, ho wever to obtain the true main sequence we used the 
model of lWrighg ( l2004bl) . This model was determined using a 
9*-order polynomial fit to the Hipparcos main sequence, de- 
scribed by: 



Mv 



,(B-y) = ^fl,(B-V)' 



(3) 



Here, My^usiB - V) is the absolute magnitude of the main 
sequence as a function of colour, a, are the coefficients describ- 
ing the polynomial, which are the following: 0.909, 6.258, - 
23.022, 125.5537, -321.1996, 485.5234, -452.3198, 249.6461, 
-73.57645 and 8.8240. As expected the majority of stars are clus- 
tered around the main sequence fit, with a few outliers above and 
below the fit. The majority of stars far below the fit are extremely 
metal-poor objects showing that metallicity can significantly al- 
ter a stars position on the HR-diagram. Conversely, high metal 
abundance will place stars above the main sequence and con- 
tribute to the spread of values around the fitted main sequence. 
The dashed lines represent the spread across the main sequence 
for a metallicity of one-third the solar metallicity (lower) and 
three times solar (upper) in a similar fashion to the analysis 
of Wright. Stars that are significantly higher than the main se- 
quence are evolving off the main sequence and moving across to 
the red giant branch. Wright has shown that the majority of stars 
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Fig. 14. The top plot shows a Hipparcos colour-magnitude dia- 
gram for all the stars in this survey except for the calibration stars 
above a B - V of 0.9. The open circles represent all the target 
stars in this colour space. The solid line is the emp irical main se- 
quence for Hipparcos stars determined in IWrigha (l2004b). with 
the dashed lines representing the spread in the main sequence 
due to metallicity (+0.3dex). The dotted line represents a differ- 
ence of 1 magnitude above the main sequence and this isolates 
evolved stars. It is a clear a number of these stars are evolving 
off the main sequence. There are also a number of a stars far be- 
low the main sequence that comprise thick-disk and halo stars. 
The filled circles represent possible Maunder Minimum stars in 
this study (upper panel). The lower plot shows the distance from 
the fitted main sequence against each stars metallicity within a 
subset of the original sample. This represents all stars where the 
surface gravity was greater than 4.0 (log(g)). The solid line rep- 
resents the best-fit straight line with the dashed lines represent- 
ing the Icr uncertainty boundaries. This fit can be used to correct 
each stars position on the main sequence for the effects of metal- 
licity. 



announced as being in a Maunder minima are actually evolved 
stars and hence the indices are not applicable to this conclu- 
sion. There are 79 stars in this study that are extremely inactive 
(log/?[jj, <-5. 10). Other authors have claimed that st ars with such 
properties are in a Maunder minimum p hase (e.g. IHenrv et all 
119961: ISaaretalJ [19981 lGravetal.ll2003l) . These stars are high- 
lighted in Fig. [14] by filled circles. It is clear that the majority 
are located high above the main sequence and are evolving onto 
the giant branch, with two below and the rest close to the main 
sequence. To determine how many of these stars are clearly on 



the main sequenc e we use the de finition of the height above the 
main sequence in lWrighg ( l2004bl) . 



{AMv)ms,b-v = Mv,Ms(B -V)-Mv 



(4) 



(AMv)ms.b-v is the difference between a stars absolute mag- 
nitude and the absolute magnitude of the main sequence fit (for 
that B - V). Eq. |4| is used to determine the height above the 
main sequence, giving an indication if a star is evolving off the 
main sequence or not. One must still be careful when working 
with solitary stars around the main sequence as their position is 
metallicity sensitive with stars of low metallicity lying below the 
main sequence and stars with high metallicities lying above the 
main sequence. Wright estimates that stars with a metallicity of 
-1-0.3 dex will lie ~0.45 magnitudes above the main sequence, 
whereas stars with a metallicity of -0.3 dex will lie ~0.45 mag- 
nitudes below the main sequence. These regions are shown as 
dashed lines in Fig.[T4|and the majority of stars in this study are 
located between these two bounds, both inactive and active stars. 
Also shown on the plot is the region that marks (AMv)ms,b-v=^ 
(dotted line) which we use as a boundary for evolutionary status. 
Stars above this line we expect to be in an evolved state of their 
evolution regardless of their metallicity. Since we have measured 
the metallicities, a correction can be made to remove the metal- 
licity component in each star's main sequence position. 

The lower plot in Figure [T4| shows the distribution of stellar 
metallicities against distance from the main sequence. A clear 
trend is seen in the plot with the metal-poor objects generally 
found below the main sequence and metal-rich objects located 
above. This sequence was cleaned by focusing on stars that show 
no significant evidence for evolution into giant status. As men- 
tioned above, all stars with a (AMv)ms,b-v >1 must be evolving 
off the main sequence and are removed from this sample. To en- 
sure this is the case and to detect any evolved objects located 
(AMv)ms,b-v <1, an estimate of the surface gravity was made 
following the technique described in Section 3.3.1. Using sur- 
face gravities we were able to remove all stars with a \og(g) of 
3.5 or 4.0, which with grid steps of 0.5, should effectively re- 
move all stars below 4.25. This is a conservative estimate but 
should select against any stars evolving off the main sequence 
and leave behind only true main sequence stars. A first fit was 
then made to this subset of stars to estimate where the main se- 
quence lies and the Icr scatter was measured and used to perform 
a secondary cut, which removed all objects outside the first scat- 
ter region. This cut was used to better constrain the position of 
the main sequence by focusing on the region that is most pop- 
ulated. The thick solid line Fig. [T4| (lower) shows the best-fit 
linear trend to this main sequence across all metallicities, after 
all cuts have been performed, accompanied by the Icr uncer- 
tainty regions (dashed lines). The trend is evident and it should 
be noted that the fit intersects the point of zero distance from the 
main sequence very close to the solar metallicity value of zero 
without any necessary fixing or adjustment. The relationship is 
described by: 



(AMv\Fe/H] = 1.050([Fe/H]) + 0.019 



(5) 



The uncertainties on the fitted coefficients in Equation. [5] are 
0.018 and 0.071 respectively, and the Icr uncertainty shown by 
the dashed lines in the figure is +0.240. Since the majority of 
stars reside within the scatter region we are confident this region 
is a good proxy for the main sequence in this sample. Simply by 
subtracting off (A Mv)[Fe/H] from Eq".|5]from the initially de- 
termined {AMv)ms,b-v in Eq".|4|we can remove the effects of 
metallicity on the spread in distance from the main sequence. 
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Table 2. Six possible Maunder Minimum stars, along with 
their Hipparcos B - V colours and V-band magnitudes, (A 
^v)ms,b-v from the main sequence after metallicity correction, 
log/?^j^ Pgj^Qg indices and metallicity values. 



HD 


B-V 


V 


AMv 


log"HK,FEROS 


[Fe/H] 


Data 
HD8446 


0.661 


8.16 


-0.15 


-5.17 


-0.35 


HD9175 


0.656 


8.01 


-0.24 


-5.12 


-0.12 


HD22104 


0.679 


8.32 


0.07 


-5.15 


0.15 


HD31827 


0.770 


8.26 


0.07 


-5.11 


-0.34 


HD84501 


0.663 


8.26 


-0.03 


-5.12 


-0.22 


HD90722 


0.724 


7.88 


-0.15 


-5.10 


-0.20 



This repositions all stars onto the main sequence irrespective 
of their metallicities and since all evolved objects have been 
removed a clean sample of possible Maunder minimum candi- 
dates can be made, based on their inactivity and distance from 
the main sequence. From the stars that are left we have selected 
those that are within the Icr boundaries of this fitted main se- 
quence. From this selection six stars remain and should rep- 
resent both main sequence stars and very inactive stars, all of 
these are listed in Table. |2] A check of any affects due to line 
strength changes between metal-poor and metal-rich objects was 
made to test if any lines within the bandpass regions were bias- 
ing the metal-rich stars towards lower activities. One absorption 
line was found inside the H bandpass region and none inside the 
K bandpass region. We interpolated this line out and recomputed 
the activity index and found that there were no real changes in 
the final logR'^^^ activities. This result was similar for the con- 
tinuum bandpass regions and therefore any line strength alter- 
ations are minimal compared with the reduction uncertainties. 
However, the metallicity does affect the intrinsic colour of stars 
and should affect the derived logR'^^ value as this is colour de- 
pendent. It is also important to note that since only one measure- 
ment has been made for these objects in this study, it is unsure 
whether or not they will remain stable over a number of years or 
will undergo any cyclical variation. If we assume a Gaussian dis- 
tribution to describe the cyclical activity variations and take the 
stars with the potential to scatter above the -5.10 activity bound- 
ary (i.e. stars within the range -5.25 to -5.10 which represents 
a maximum cycle of 0.15 in logR[^^ and is consistent with the 
Sun's cycle) then four of the six stars are expected to be nor- 
mal stars undergoing their usual activity cycle at the Icr level 
of confidence. Of the entire sample of 353 stars, this leads to a 
fraction of ~2+l% and hence the Sun should spend only a few 
percent of its main sequence lifetime in a Maunder minimum 
like state. It should be noted that the sample is biased towards 
the most metal-rich objects in the solar neighbourhood. This se- 
lection will skew the final estimate of Maunder minimum stars 
as it should be biased towards younger more active objects. 

5.4. Benchmark Target List 

The main aim of this study was to determine accurate chromo- 
spheric activity and metallicity values for a number of close-by, 
bright stars in the southern hemisphere. The metal-rich subset of 
these will provide excellent candidates for planet search studies 
searching for hot giant planets. Figure [12] (upper) shows the dis- 
tribution of metallicities in this sample. About 30% of stars have 
metallicities greater than 0.1 dex, of which 15% have metal- 
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Fig. 15. Iron abundance as a function of chromospheric activity 
for the stars in this sample. The filled circles represent all stars 
that will compose our metal-rich target list. While most stars 
above -hO. 1 dex are inactive, a few have activities above -4.5. 



licities greater than 0.2 dex and it is these stars that represent 
the best stars for radial-velocity follow-up. First epoch radial- 
velocities have already been acquired for a number of these stars 
in ESO Period 76 with the HARPS instrument. 

The chromospheric activities also play an important role in 
the nature of the objects that can be detected. Figure [15] shows 
the distribution of activities across all metallicity space. Only a 
few of the metal-rich candidates are active and will be removed 
from the final candidate list, whereas most of the stars above 
-t-O.ldex in [Fe/H] have logR'^^ values less than -4.5. Since the 
level of radial-velocity noise (aka. jitter) is coiTelated with ac- 
tivity (Saaretal. 1998; Santos et al. 2000; Wright 2005) these 
inactive stars allow higher precision radial-velocities to be effi- 
ciently measured. The potential to access lower amplitude sig- 
nals allows access to objects of lower masses and as the mass 
function rises towards lower masses following a power law of 
-1.8 (Grether & Lineweaver 2006), the probability of planet de- 
tection significantly increases. The final target list will be a sub- 
set of stars taken from Table. [4] following the criteria explained 
above. 



6. Conclusions 

Using the high-resolution echelle data from the FEROS instru- 
ment we extracted chromospheric activities and metallicities for 
353 bright stars in the Southern hemisphere. We have shown the 
activity distribution can be described well by a bimodal function 
which agrees well with previous work in both the northern and 
southern hemispheres for such stars. The sample is dominated by 
G stars (rEFF~6000K) and the activity distribution peaks around 
-5.0. The K stars have an active to inactive peak ratio close to 1, 
compared with the G stars peak-to-peak ratio of 0.2. The metal- 
licity distribution peaks around 0.1 dex, contrary to unbiased 
surveys which peak around -0.1 dex. 116 stars have metallici- 
ties greater than 0. 1 dex and are not found on other planet search 
lists. It is these stars that constitute our planned metal-rich planet 
search project. The uncertainties follow a Gaussian distribution 
peaking at ±0.085 dex with a high uncertainty tail of objects due 
to stars with low S/N spectra. The val ues are in agreement wit h 
similar works in the literature, such as lValenti & Fischen (l2005h . 
with an RMS scatter of -0.05 dex. 
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16 J.S. Jenkins et al.: Metallicides & Activities of Southern Stars 

Table 3. Derived activity and metallicity values for all stars used to calibrate onto the Mt. Wilson system of measurements. HD 
Number refers to the Henry Draper catalogue identifier. Two stars have no HD catalogue reference and are catalogued by their 
Hipparcos identifier (HIP). Johnson B-V, V taken from the Hipparcos catalogue. S'FEROS.a is before calibration onto the Mt. Wilson 
and ^FEROS.b is after calibra tion. 5mw is the Mt. Wilson index used to perform t his calibration, take n from column 5 of table 1 and 
columns 4-7 of table 3 from lDuncan et al] (Il99lh . logTJ^j, ^^ are derived from Duncan et al.l (Il99ll) . Note that the photon-counting 
uncertainty in the derived activities are indicative only; based on measurements of t Ceti a random uncertainty of t 3% must be 
taken into account in all instances. [Fe/H] is used to denote iron abundances (herein referred to as metallicity). 



HD 


B-V 


V 


^'pEROS.a 


■SpEROS.b 


5mw 


'°8"hk,feros 


'°8"hk,mw 


[Fe/H] 


Calibrators 


















HD1835 


0.659 


6.39 


0.291 


0.350±0.003 


0.347 


-4.43 


-4.44 


0.21+0.03 


HD3443AB 


0.715 


5.57 


0.142 


0.180±0.001 


0.186 


-4.92 


-4.89 


-0.01+0.05 


HD3795 


0.718 


6.14 


0.119 


0.153±0.001 


0.156 


-5.06 


-5.04 


-0.65+0.05 


HD9562 


0.639 


5.75 


0.106 


0.138±0.001 


0.144 


-5.16 


-5.11 


0.19+0.14 


HD 10700 


0.727 


3.49 


0.140 


0.178±0.001 


0.173 


-4.93 


-4.95 


-0.38+0.04 


HD 10700 


0.727 


3.49 


0.142 


0.180±0.001 


0.173 


-4.92 


-4.95 


-0.35+0.04 


HD22484 


0.575 


4.29 


0.118 


0.153±0.001 


0.147 


-5.02 


-5.06 


-0.14+0.02 


HD30495 


0.632 


5.49 


0.234 


0.286+0.002 


0.292 


-4.54 


-4.54 


0.03+0.03 
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Table 4. S'feros and log^J^j^ pgj^Qs for bright solar-type stars in the southern hemisphere. HD Number refers to the Henry Draper 
catalogue identifier. Johnson B-V, V and Spec Type are all taken from the Hipparcos catalogue. 5feros is the FEROS activity 
index after calibration onto the Mt. Wilson system of measurements. log^HKFEROS ^^^ ^^^^ activity values generated following the 
iNoves et al] (Il984al) methodology. [Fe/H] is the iron abundance. 



HD/HIP 


B-V 


V 


■Sferos 


log"HK.FEROS 


[Fe/H] 


Measurements 












HD2595 


0.840 


9.34 


0.174+0.001 


-5.00 


-0.40+0.07 


HD3569 


0.846 


9.21 


0.176+0.001 


-4.99 


-0.29+0.08 


HD5185 


0.722 


8.92 


0.239+0.002 


-4.71 


-0.08+0.07 


HD6348 


0.801 


9.15 


0.177+0.001 


-4.97 


-0.50+0.06 


HD6558 


0.606 


8.20 


0.141+0.001 


-5.12 


0.25+0.06 


HD6790 


0.561 


9.37 


0.168+0.001 


-4.92 


0.20+0.03 


HD7112 


0.672 


8.80 


0.382+0.003 


-4.39 


-0.04+0.04 


HD7320 


0.685 


8.95 


0.163+0.001 


-4.99 


-0.28+0.06 


HD7621 


0.762 


9.03 


0.142+0.001 


-5.13 


-0.01+0.07 


HD7950 


0.797 


8.72 


0.139+0.001 


-5.16 


0.17+0.07 


HD8038 


0.701 


8.41 


0.210+0.002 


-4.79 


0.22+0.06 


HD8389 


0.900 


7.85 


0.152+0.001 


-5.11 


0.40+0.08 


HD8446 


0.661 


8.16 


0.137+0.001 


-5.17 


0.28+0.04 


HD8765 


0.707 


8.14 


0.180+0.001 


-4.91 


0.23+0.06 


HD9279 


0.630 


8.15 


0.164+0.001 


-4.97 


-0.08+0.04 


HD9796 


0.780 


8.77 


0.179+0.001 


-4.95 


-0.37+0.07 


HD 10008 


0.797 


7.66 


0.468+0.004 


-4.39 


-0.02+0.06 


HD10188 


0.676 


8.94 


0.140+0.001 


-5.14 


0.30+0.06 


HD10278 


0.695 


9.13 


0.279+0.002 


-4.59 


0.15+0.07 


HD11352 


0.644 


8.61 


0.308+0.003 


-4.50 


-0.12+0.04 


HD 11683 


0.887 


9.15 


0.508+0.004 


-4.45 


-0.26+0.08 


HD 12754 


0.566 


8.66 


0.206+0.002 


-4.74 


-0.11+0.14 


HD12873 


0.795 


8.95 


0.151+0.001 


-5.09 


-0.14+0.07 


HD12889 


0.885 


8.55 


0.167+0.001 


-5.05 


-0.37+0.07 


HD13147 


0.695 


9.21 


0.136+0.001 


-5.18 


0.16+0.07 


HD13350 


0.766 


8.94 


0.131+0.001 


-5.21 


0.34+0.08 


HD 13977 


0.880 


9.11 


0.244+0.002 


-4.82 


0.07+0.10 


HD 14655 


0.684 


8.13 


0.138+0.001 


-5.16 


-0.01+0.05 


HD15191 


0.704 


8.65 


0.193+0.002 


-4.85 


0.07+0.06 


HD15507 


0.670 


8.66 


0.178+0.001 


-4.90 


0.22+0.06 


HD17194 


0.803 


8.71 


0.122+0.001 


-5.28 


0.06+0.08 


HD 18241 


0.672 


8.66 


0.346+0.003 


-4.45 


-0.10+0.03 


HD 18599 


0.876 


8.99 


0.527+0.004 


-4.42 


-0.01+0.12 


HD 18599 


0.876 


8.99 


0.544+0.005 


-4.40 


-0.01+0.13 


HD 18708 


0.651 


8.42 


0.136+0.001 


-5.18 


0.18+0.05 


HD 18754 


0.736 


8.44 


0.131+0.001 


-5.22 


0.16+0.06 


HD 19493 


0.695 


8.14 


0.140+0.001 


-5.15 


0.23+0.06 


HD19773 


0.551 


8.06 


0.177+0.001 


-4.86 


0.22+0.03 


HD20605 


0.766 


9.22 


0.176+0.001 


-4.95 


-0.49+0.06 


HD22104 


0.679 


8.32 


0.139+0.001 


-5.15 


0.31+0.05 


HD22177 


0.717 


7.78 


0.152+0.001 


-5.06 


0.16+0.06 


HD22282 


0.789 


8.54 


0.164+0.001 


-5.02 


0.16+0.07 


HD22472 


0.784 


8.85 


0.229+0.002 


-4.78 


0.04+0.09 


HD23065 


0.751 


8.25 


0.156+0.001 


-5.05 


-0.59+0.04 


HD23398 


0.788 


8.43 


0.133+0.001 


-5.19 


0.44+0.07 


HD24091 


0.817 


8.89 


0.485+0.004 


-4.39 


-0.06+0.03 


HD25015 


0.899 


8.87 


0.537+0.005 


-4.44 


0.04+0.07 


HD25061 


0.843 


9.27 


0.238+0.002 


-4.80 


0.12+0.08 


HD26071 


0.742 


8.89 


0.157+0.001 


-5.04 


0.21+0.06 


HD26965 


0.820 


4.43 


0.181+0.001 


-4.96 


-0.27+0.09 


HD27905 


0.627 


7.82 


0.174+0.001 


-4.91 


-0.08+0.06 


HD28013 


0.701 


8.38 


0.186+0.001 


-4.88 


0.01+0.05 


HD28255 


0.659 


6.28 


0.171+0.001 


-4.94 


0.08+0.04 


HD28255 


0.659 


6.28 


0.179+0.001 


-4.89 


0.05+0.04 


HD29231 


0.776 


7.61 


0.246+0.002 


-4.73 


0.11+0.07 


HD29813 


0.621 


7.74 


0.173+0.001 


-4.91 


-0.05+0.05 


HD30501 


0.875 


7.58 


0.334+0.003 


-4.65 


-0.06+0.08 


HD31827 


0.770 


8.26 


0.147+0.001 


-5.11 


0.39+0.07 


HD33449 


0.681 


8.49 


0.180+0.001 


-4.90 


-0.47+0.11 


HD33725 


0.799 


8.04 


0.220+0.002 


-4.82 


-0.15+0.07 


HD34688 


0.834 


9.08 


0.203+0.002 


-4.89 


-0.21+0.07 


HD34962 


0.733 


8.27 


0.383+0.003 


-4.44 


-0.04+0.05 


HD36169 


0.702 


8.32 


0.170+0.001 


-4.96 


-0.09+0.05 
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HD37634 


0.647 


8.91 


0.186±0.001 


-4.86 


-0.51±0.03 


HD38459 


0.861 


8.52 


0.471 ±0.004 


-4.46 


0.13±0.06 


HD38467 


0.672 


8.27 


0.147±0.001 


-5.09 


0.22±0.05 


HD39833 


0.629 


7.65 


0.203±0.002 


-4.78 


0.15±0.04 


HD39997 


0.559 


8.48 


0.161±0.001 


-4.96 


0.06±0.03 


HD40293 


0.720 


9.10 


0.206±0.002 


-4.81 


0.13±0.06 


HD42538 


0.593 


8.21 


0.140±0.001 


-5.13 


0.22±0.04 


HD42719 


0.657 


7.51 


0.138±0.001 


-5.16 


0.24±0.04 


HD42936 


0.900 


9.10 


0.157±0.001 


-5.09 


0.18±0.09 


HD43197 


0.817 


8.98 


0.146±0.001 


-5.12 


0.37±0.09 


HD43306 


0.835 


8.70 


0.168±0.001 


-5.02 


0.10±0.08 


HD44310 


0.837 


8.68 


0.187±0.001 


-4.94 


-0.02±0.08 


HD44569 


0.619 


7.73 


0.176±0.001 


-4.90 


-0.21±0.03 


HD45133 


0.741 


8.41 


0.147±0.001 


-5.10 


0.30±0.06 


HD46894 


0.778 


7.92 


0.143±0.001 


-5.13 


0.02±0.06 


HD47186 


0.714 


7.63 


0.163±0.001 


-5.00 


0.21±0.06 


HD47252 


0.781 


8.27 


0.192±0.002 


-4.90 


-0.45±0.06 


HD48265 


0.747 


8.05 


0.128±0.001 


-5.24 


0.29±0.05 


HD48611 


0.742 


8.71 


0.228±0.002 


-4.76 


-0.32±0.05 


HD49035 


0.755 


8.53 


0.225±0.002 


-4.77 


0.29±0.06 


HD49866 


0.663 


8.03 


0.144±0.001 


-5.11 


0.10±0.09 


HD50652 


0.691 


8.65 


0.152±0.001 


-5.06 


0.24±0.06 


HD52756 


0.899 


8.47 


0.246±0.002 


-4.83 


0.09±0.09 


HD55524 


0.722 


9.42 


0.148±0.001 


-5.09 


0.24±0.06 


HD56259 


0.760 


8.78 


0.145±0.001 


-5.12 


0.14±0.07 


HD56413 


0.754 


8.91 


0.169±0.001 


-4.98 


0.24±0.06 


HD56662 


0.604 


7.67 


0.164±0.001 


-4.96 


-0.16±0.03 


HD56957 


0.701 


7.57 


0.138±0.001 


-5.16 


0.26±0.05 


HD58111 


0.875 


8.83 


0.211 ±0.002 


-4.90 


-0.04±0.08 


HD59100 


0.634 


8.17 


0.172±0.001 


-4.92 


-0.15±0.04 


HD61033 


0.724 


7.59 


0.401±0.003 


-4.40 


-0.32±0.06 


HD61475 


0.853 


8.79 


0.380±0.003 


-4.56 


0.22±0.07 


HD62128 


0.717 


7.53 


0.146±0.001 


-5.10 


0.30±0.07 


HD62549 


0.612 


7.72 


0.168±0.001 


-4.94 


0.04±0.06 


HD63765 


0.745 


8.10 


0.226±0.002 


-4.77 


-0.13±0.06 


HD64010 


0.815 


9.26 


0.418±0.004 


-4.47 


-0.19±0.05 


HD65562 


0.869 


8.81 


0.213±0.002 


-4.89 


-0.17±0.08 


HD66221 


0.721 


8.06 


0.161±0.001 


-5.01 


0.16±0.06 


HD66340 


0.831 


9.25 


0.201±0.002 


-4.89 


0.04±0.07 


HD66428 


0.715 


8.25 


0.152±0.001 


-5.06 


0.25±0.05 


HD66653 


0.655 


7.52 


0.210±0.002 


-4.76 


0.15±0.06 


HD67200 


0.595 


7.71 


0.151±0.001 


-5.04 


0.30±0.04 


HD68475 


0.887 


8.78 


0.324±0.003 


-4.67 


-0.07±0.09 


HD68785 


0.619 


8.19 


0.196±0.002 


-4.80 


-0.21±0.03 


HD68916 


0.820 


9.39 


0.211 ±0.002 


-4.86 


-0.34±0.06 


HD69611 


0.584 


7.74 


0.165±0.001 


-4.94 


-0.61±0.07 


HD69721 


0.826 


8.80 


0.167±0.001 


-5.02 


0.15±0.09 


HD70081 


0.660 


8.92 


0.147±0.001 


-5.09 


0.16±0.05 


HD71251 


0.851 


9.19 


0.371±0.003 


-4.57 


0.01±0.07 


HD73536 


0.698 


8.19 


0.164±0.001 


-4.99 


-0.06±0.05 


HD75288 


0.673 


8.51 


0.159±0.001 


-5.01 


0.09±0.05 


HD75519 


0.651 


7.83 


0.401±0.003 


-4.35 


0.13±0.05 


HD76849 


0.843 


9.06 


0.164±0.001 


-5.04 


0.05±0.07 


HD77134 


0.537 


7.80 


0.179±0.001 


-4.84 


-0.28±0.01 


HD78130 


0.674 


8.66 


0.283±0.002 


-4.57 


0.17±0.04 


HD78286 


0.682 


8.59 


0.148±0.001 


-5.09 


0.25±0.05 


HD78663 


0.767 


8.62 


0.184±0.001 


-4.92 


-0.46±0.06 


HD78747 


0.575 


7.72 


0.170±0.001 


-4.91 


-0.69±0.01 


HD81044 


0.801 


8.85 


0.377±0.003 


-4.51 


-0.41±0.07 


HD82977 


0.790 


9.08 


0.305±0.003 


-4.61 


-0.05±0.05 


HD83517 


0.664 


7.74 


0.172±0.001 


-4.93 


-0.23±0.06 


HD84501 


0.663 


8.26 


0.144±0.001 


-5.12 


0.09±0.04 


HD84742 


0.558 


7.93 


0.168±0.001 


-4.91 


-0.16±0.43 


HD85228 


0.894 


7.91 


0.230±0.002 


-4.86 


-0.03±0.09 


HD85404 


0.591 


7.78 


0.152±0.001 


-5.03 


0.01±0.07 


HD86006 


0.708 


8.17 


0.155±0.001 


-5.04 


0.21±0.05 


HD86171 


0.746 


8.53 


0.260±0.002 


-4.67 


-0.17±0.06 


HD87007 


0.840 


8.82 


0.156±0.001 


-5.07 


0.16±0.08 


HD87978 


0.687 


8.15 


0.362±0.003 


-4.43 


0.02±0.11 


HD88261 


0.590 


8.07 


0.174±0.001 


-4.89 


-0.47±0.21 


HD88864 


0.637 


7.89 


0.172±0.001 


-4.92 


-0.00±0.07 


HD90028 


0.656 


8.32 


0.150±0.001 


-5.07 


0.16±0.04 
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HD90520 


0.643 


7.51 


0.161±0.001 


-4.99 


0.17+0.14 


HD90722 


0.724 


7.88 


0.146±0.001 


-5.10 


0.28+0.06 


HD90812 


0.813 


8.86 


0.182±0.001 


-4.95 


-0.33+0.06 


HD91682 


0.696 


8.51 


0.153±0.001 


-5.06 


0.16+0.05 


HD92213 


0.802 


9.41 


0.192±0.002 


-4.91 


-0.43+0.07 


HD93528 


0.831 


8.39 


0.514±0.004 


-4.38 


-0.11+0.03 


HD93626 


0.536 


7.88 


0.177±0.001 


-4.86 


0.06+0.02 


HD93849 


0.559 


7.86 


0.147±0.001 


-5.06 


0.27+0.03 


HD93932 


0.615 


7.53 


0.154±0.001 


-5.02 


0.10+0.05 


HD94151 


0.718 


7.84 


0.184±0.001 


-4.90 


0.05+0.05 


HD94268 


0.557 


7.63 


0.175±0.001 


-4.88 


-0.39+0.02 


HD94270 


0.579 


7.76 


0.163±0.001 


-4.95 


0.00+0.03 


HD94482 


0.562 


8.13 


0.157±0.001 


-4.99 


-0.12+0.01 


HD94690 


0.710 


8.22 


0.158±0.001 


-5.02 


0.11+0.05 


HD95136 


0.660 


8.65 


0.157±0.001 


-5.02 


0.17+0.04 


HD95338 


0.878 


8.62 


0.185±0.001 


-4.98 


0.06+0.08 


HD96020 


0.710 


8.80 


0.151±0.001 


-5.07 


0.20+0.06 


HD96494 


0.809 


8.87 


0.344±0.003 


-4.57 


0.14+0.06 


HD96941 


0.736 


8.69 


0.189±0.002 


-4.89 


-0.25+0.05 


HD97038 


0.661 


7.70 


0.155±0.001 


-5.03 


0.15+0.05 


HD97356 


0.637 


8.15 


0.166±0.001 


-4.95 


0.01+0.04 


HD97517 


0.613 


8.11 


0.165±0.001 


-4.95 


-0.35+0.02 


HD98356 


0.828 


8.73 


0.228±0.002 


-4.81 


0.07+0.07 


HD98640 


0.765 


8.53 


0.335±0.003 


-4.54 


-0.02+0.03 


HD98649 


0.658 


8.00 


0.169±0.001 


-4.95 


-0.06+0.04 


HD99565 


0.748 


7.65 


0.254±0.002 


-4.69 


-0.37+0.01 


HD100555 


0.726 


8.17 


0.207±0.002 


-4.81 


0.07+0.06 


HD 100922 


0.804 


8.86 


0.199±0.002 


-4.88 


-0.42+0.04 


HD101093 


0.561 


7.61 


0.185±0.001 


-4.83 


-0.12+0.35 


HD101197 


0.659 


8.73 


0.171±0.001 


-4.93 


0.19+0.06 


HD101348 


0.711 


7.82 


0.138±0.001 


-5.17 


0.23+0.06 


HD102196 


0.569 


8.15 


0.146±0.001 


-5.07 


0.13+0.03 


HD 102346 


0.551 


8.15 


0.322±0.003 


-4.42 


0.25+0.05 


HD102361 


0.556 


8.16 


0.157±0.001 


-4.98 


0.12+0.03 


HD103197 


0.860 


9.40 


0.158±0.001 


-5.07 


-0.12+0.04 


HD103673 


0.705 


8.59 


0.359±0.003 


-4.45 


0.02+0.06 


HD 104006 


0.780 


8.91 


0.179±0.001 


-4.95 


-0.61+0.05 


HD 104263 


0.753 


8.23 


0.169±0.001 


-4.98 


-0.04+0.06 


HD104576 


0.708 


8.53 


0.422±0.004 


-4.36 


-0.55+0.05 


HD104988 


0.754 


8.16 


0.166±0.001 


-5.00 


-0.24+0.05 


HD 105750 


0.739 


8.96 


0.157±0.001 


-5.04 


0.20+0.03 


HD105837 


0.570 


7.52 


0.185±0.001 


-4.83 


-0.34+0.04 


HD 105901 


0.626 


8.20 


0.168±0.001 


-4.94 


0.03+0.05 


HD106275 


0.892 


8.63 


0.207±0.002 


-4.92 


-0.14+0.02 


HD 106937 


0.762 


8.50 


0.137±0.001 


-5.17 


0.18+0.07 


HD107181 


0.737 


8.43 


0.140±0.001 


-5.15 


0.26+0.07 


HD107859 


0.840 


9.42 


0.182±0.001 


-4.96 


-0.39+0.05 


HD108341 


0.920 


9.38 


0.181±0.001 


-5.02 


-0.12+0.05 


HD 108446 


0.872 


8.72 


0.331±0.003 


-4.65 


-0.19+0.06 


HD108523 


0.710 


8.31 


0.164±0.001 


-4.99 


-0.00+0.06 


HD108682 


0.784 


9.06 


0.216±0.002 


-4.82 


-0.43+0.04 


HD108953 


0.800 


9.26 


0.174±0.001 


-4.98 


0.29+0.07 


HD 109423 


0.898 


8.82 


0.423±0.004 


-4.55 


0.06+0.01 


HD109570 


0.692 


8.04 


0.217±0.002 


-4.76 


-0.05+0.05 


HD 109908 


0.583 


8.34 


0.294+0.002 


-4.49 


0.13+0.05 


HD109930 


0.894 


8.52 


0.173±0.001 


-5.03 


0.41+0.07 


HDl 10255 


0.882 


9.34 


0.191±0.002 


-4.96 


-0.14+0.04 


HD 11 0605 


0.730 


9.13 


0.158±0.001 


-5.03 


0.23+0.07 


HDl 10652 


0.819 


9.47 


0.171±0.001 


-4.99 


-0.39+0.05 


HDl 10979 


0.654 


8.05 


0.163±0.001 


-4.98 


0.09+0.06 


HDl 11096 


0.703 


7.69 


0.153±0.001 


-5.05 


0.25+0.06 


HD111431 


0.627 


8.01 


0.144±0.001 


-5.11 


0.04+0.05 


HDl 11777 


0.616 


8.48 


0.173±0.001 


-4.91 


-0.44+0.03 


HDl 12037 


0.658 


8.53 


0.185±0.001 


-4.87 


-0.11+0.05 


HD112121 


0.728 


8.93 


0.149±0.001 


-5.09 


0.13+0.06 


HDl 12758 


0.769 


7.54 


0.176±0.001 


-4.96 


-0.50+0.04 


HDl 12863 


0.779 


8.70 


0.407±0.003 


-4.45 


-0.30+0.03 


HDl 13933 


0.746 


8.75 


0.158±0.001 


-5.03 


0.01+0.06 


HDl 17579 


0.721 


8.87 


0.169±0.001 


-4.97 


-0.14+0.05 


HDl 19586 


0.870 


9.12 


0.171±0.001 


-5.02 


-0.75+0.01 


HDl 19607 


0.876 


9.21 


0.409±0.003 


-4.54 


-0.16+0.05 
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HD 120067 


0.797 


8.70 


0.328±0.003 


-4.58 


-0.21+0.05 


HD 123433 


0.875 


9.25 


0.182±0.001 


-4.99 


-0.17+0.06 


HD 124595 


0.610 


8.76 


0.152±0.001 


-5.04 


0.06+0.05 


HD126535 


0.840 


8.86 


0.479±0.004 


-4.43 


0.13+0.07 


HD 126828 


0.533 


8.25 


0.142±0.001 


-5.10 


-0.06+0.04 


HD 127423 


0.569 


8.53 


0.239±0.002 


-4.63 


0.13+0.05 


HD 128356 


0.685 


8.29 


0.218+0.002 


-4.75 


-0.78+0.03 


HD 128985 


0.833 


9.24 


0.172±0.001 


-5.00 


-0.37+0.04 


HD 129445 


0.756 


8.80 


0.175±0.001 


-4.95 


0.27+0.08 


HD143120 


0.780 


7.52 


0.125±0.001 


-5.26 


0.23+0.07 


HD143361 


0.773 


9.20 


0.136±0.001 


-5.18 


0.06+0.07 


HD144550 


0.682 


8.61 


0.137±0.001 


-5.17 


0.22+0.06 


HD 144848 


0.648 


8.66 


0.140±0.001 


-5.15 


0.28+0.07 


HD 144899 


0.660 


8.97 


0.132±0.001 


-5.22 


0.36+0.20 


HD146233 


0.652 


5.49 


0.167+0.001 


-4.96 


0.30+0.03 


HD 146434 


0.700 


8.14 


0.154+0.001 


-5.05 


0.09+0.06 


HD146856 


0.870 


9.40 


0.332+0.003 


-4.64 


-0.23+0.06 


HD146856 


0.870 


9.40 


0.357+0.003 


-4.61 


-0.23+0.07 


HD 147092 


0.798 


9.41 


0.127+0.001 


-5.23 


-0.32+0.02 


HD148156 


0.560 


7.69 


0.145+0.001 


-5.08 


0.22+0.15 


HD148319 


0.645 


8.58 


0.284+0.002 


-4.55 


0.28+0.06 


HD148577 


0.664 


7.93 


0.137+0.001 


-5.17 


0.06+0.11 


HD148577 


0.664 


7.93 


0.139+0.001 


-5.15 


0.08+0.32 


HD149189 


0.653 


8.56 


0.133+0.001 


-5.21 


0.27+0.16 


HD149194 


0.559 


8.68 


0.153+0.001 


-5.01 


0.13+0.04 


HD149616 


0.561 


8.21 


0.152+0.001 


-5.02 


0.08+0.01 


HD 149724 


0.783 


7.85 


0.151+0.001 


-5.08 


0.30+0.07 


HD149782 


0.759 


9.16 


0.188+0.002 


-4.90 


0.28+0.11 


HD149782 


0.759 


9.16 


0.189+0.002 


-4.90 


0.25+0.10 


HD 150437 


0.683 


7.84 


0.127+0.001 


-5.27 


0.39+0.09 


HD 150936 


0.740 


8.76 


0.172+0.001 


-4.96 


0.16+0.07 


HD 152079 


0.711 


9.18 


0.152+0.001 


-5.06 


0.32+0.09 


HD 154221 


0.640 


8.60 


0.157+0.001 


-5.01 


0.23+0.06 


HD 154672 


0.713 


8.21 


0.138+0.001 


-5.17 


0.21+0.07 


HD 157691 


0.642 


8.37 


0.161+0.001 


-4.99 


-0.30+0.05 


HD 157798 


0.684 


8.15 


0.135+0.001 


-5.19 


0.00+0.11 


HD 15 8469 


0.556 


7.93 


0.136+0.001 


-5.16 


0.11+0.04 


HD162816 


0.630 


8.48 


0.139+0.001 


-5.15 


0.08+0.07 


HD165155 


1.018 


9.36 


0.151+0.001 


-5.20 


0.10+0.08 


HD 165204 


0.770 


9.14 


0.134+0.001 


-5.19 


0.30+0.07 


HD 165920 


0.842 


7.91 


0.151+0.001 


-5.09 


0.19+0.08 


HD166101 


0.793 


8.83 


0.410+0.003 


-4.46 


-0.22+0.06 


HD 166745 


0.763 


8.49 


0.137+0.001 


-5.17 


0.23+0.08 


HD 170706 


0.789 


9.48 


0.134+0.001 


-5.19 


0.24+0.09 


HD173859 


0.737 


8.82 


0.372+0.003 


-4.46 


0.11+0.06 


HD173872 


0.889 


8.45 


0.180+0.001 


-5.00 


-0.54+0.03 


HD 174541 


0.631 


8.59 


0.143+0.001 


-5.11 


-0.13+0.05 


HD175129 


0.622 


8.47 


0.308+0.003 


-4.48 


-0.17+0.10 


HD175167 


0.751 


8.01 


0.134+0.001 


-5.19 


-0.17+0.14 


HD176157 


0.829 


8.39 


0.204+0.002 


-4.88 


-0.20+0.06 


HD 178340 


0.757 


8.17 


0.142+0.001 


-5.13 


0.24+0.08 


HD178340 


0.757 


8.17 


0.144+0.001 


-5.12 


0.32+0.07 


HD178787 


0.866 


9.16 


0.261+0.002 


-4.77 


-0.12+0.05 


HD185181 


0.872 


9.11 


0.891+0.008 


-4.17 


0.09+0.09 


HD 185679 


0.686 


8.72 


0.160+0.001 


-5.00 


0.14+0.07 


HD186194 


0.697 


8.64 


0.145+0.001 


-5.11 


0.18+0.07 


HD 186265 


0.790 


9.26 


0.139+0.001 


-5.15 


0.32+0.09 


HD 187620 


0.692 


8.43 


0.147+0.001 


-5.09 


-0.03+0.06 


HD 187694 


0.746 


9.36 


0.161+0.001 


-5.02 


-0.21+0.08 


HD 188298 


0.657 


8.46 


0.240+0.002 


-4.67 


0.13+0.06 


HD188581 


0.614 


8.35 


0.307+0.003 


-4.48 


-0.21+0.04 


HD 188903 


0.565 


8.27 


0.142+0.001 


-5.10 


-0.51+0.63 


HD 189627 


0.556 


8.11 


0.128+0.001 


-5.26 


0.37+0.04 


HD190125 


0.708 


8.93 


0.227+0.002 


-4.74 


0.22+0.09 


HD190613 


0.629 


8.12 


0.139+0.001 


-5.15 


0.13+0.06 


HD 190647 


0.743 


7.78 


0.132+0.001 


-5.21 


0.23+0.07 


HD191122 


0.643 


8.65 


0.147+0.001 


-5.09 


0.37+0.09 


HD191760 


0.668 


8.26 


0.137+0.001 


-5.17 


0.29+0.07 


HD 192771 


0.723 


8.64 


0.160+0.001 


-5.02 


-0.15+0.05 


HD 193690 


0.790 


9.37 


0.185+0.001 


-4.93 


0.24+0.08 


HD 193 844 


0.888 


9.24 


0.238+0.002 


-4.84 


-0.20+0.01 


HD 193995 


0.716 


8.51 


0.132+0.001 


-5.21 


0.23+0.07 
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HD 194490 


0.637 


8.99 


0.137±0.001 


-5.17 


0.28+0.08 


HD 195 145 


0.738 


8.66 


0.161±0.001 


-5.02 


0.26+0.08 


HD 197499 


0.576 


8.62 


0.149±0.001 


-5.05 


0.33+0.07 


HD200869 


0.851 


9.38 


0.140±0.001 


-5.15 


0.16+0.09 


HD201378 


0.625 


7.62 


0.302±0.003 


-4.50 


-0.14+0.04 


HD201757 


0.714 


8.19 


0.134±0.001 


-5.20 


0.15+0.06 


HD203384 


0.762 


8.03 


0.163±0.001 


-5.01 


0.22+0.07 


HD204807 


0.717 


8.28 


0.163±0.001 


-5.00 


-0.20+0.06 


HD204941 


0.878 


8.45 


0.188±0.002 


-4.97 


-0.25+0.06 


HD205855 


0.844 


8.62 


0.169±0.001 


-5.02 


-0.17+0.07 


HD205871 


0.641 


8.39 


0.158±0.001 


-5.01 


-0.25+0.05 


HD206116 


0.576 


7.61 


0.133±0.001 


-5.20 


0.25+0.33 


HD206172 


0.672 


8.54 


0.177±0.001 


-4.91 


-0.08+0.06 


HD206683 


0.657 


8.32 


0.148±0.001 


-5.08 


0.33+0.06 


HD207832 


0.691 


8.78 


0.207±0.002 


-4.80 


0.16+0.07 


HD208672 


0.641 


8.30 


0.216±0.002 


-4.73 


0.18+0.07 


HD209713 


0.754 


9.07 


0.289±0.002 


-4.62 


0.03+0.08 


HD2 12521 


0.770 


8.63 


0.419+0.004 


-4.42 


0.02+0.06 


HD2 16777 


0.659 


8.01 


0.166±0.001 


-4.96 


-0.21+0.05 


HD2 18249 


0.879 


9.33 


0.196±0.002 


-4.94 


-0.73+0.01 


HD2 18249 


0.879 


9.33 


0.203±0.002 


-4.92 


-0.72+0.01 


HD2 18750 


0.880 


9.25 


0.163±0.001 


-5.06 


-0.17+0.07 


HD2 18960 


0.667 


9.08 


0.138±0.001 


-5.16 


0.21+0.07 


HD219011 


0.727 


8.90 


0.135±0.001 


-5.19 


0.26+0.08 


HD219556 


0.759 


9.08 


0.145±0.001 


-5.11 


0.04+0.08 


HD219556 


0.759 


9.08 


0.145±0.001 


-5.11 


0.05+0.08 


HD220256 


0.853 


8.56 


0.166±0.001 


-5.03 


-0.25+0.06 


HD220981 


0.740 


8.75 


0.142+0.001 


-5.13 


0.18+0.07 


HD221343 


0.657 


8.37 


0.232+0.002 


-4.69 


0.16+0.07 


HD221575 


0.890 


8.81 


0.523+0.004 


-4.44 


0.28+0.03 


HD221954 


0.750 


8.72 


0.139+0.001 


-5.16 


0.27+0.07 


HD222013 


0.809 


9.22 


0.161+0.001 


-5.04 


0.05+0.08 


HD222342 


0.654 


8.87 


0.180+0.001 


-4.89 


0.11+0.07 


HD222910 


0.765 


9.24 


0.132+0.001 


-5.21 


0.12+0.09 


HD223078 


0.882 


8.93 


0.210+0.002 


-4.91 


-0.73+0.01 


HD223315 


0.742 


8.78 


0.208+0.002 


-4.82 


0.25+0.08 


HD223713 


0.794 


9.43 


0.162+0.001 


-5.03 


-0.43+0.05 


HD224065 


0.680 


8.86 


0.180+0.001 


-4.90 


-0.18+0.07 


HD224433 


0.748 


8.84 


0.165+0.001 


-5.00 


0.12+0.07 


HD224538 


0.581 


8.06 


0.148+0.001 


-5.06 


0.39+0.06 


HD224828 


0.642 


8.57 


0.168+0.001 


-4.95 


-0.23+0.07 


HD224908 


0.620 


8.64 


0.186+0.001 


-4.84 


0.00+0.06 


HD302554 


0.705 


8.83 


0.332+0.003 


-4.50 


-0.02+0.06 


HIP37727 


0.700 


7.55 


0.288+0.003 


-4.58 


-0.57+0.06 



